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en liaison avec leur structure 
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SoMMAIRE. — On expose les résultats fournis par Vélude expérimentale d’un certain nombre de propriétés optiques des couches minces 


transparentes déposées par évaporation thermique, et on montre comment on peut utiliser ces résultats pour délerminer la structure 

de ces couches. Les mesures effectuées ont surtout été des mesures spectrophotométriques en incidence normale et en incidence 

ee ; on a éludié de facon approfondie les phénoménes qui se produisent pour des incidences situées dans la zone de réflexion 
e. 

Tous les résultats oblenus concourent a prouver que les couches minces étudiées ont une structure beaucoup plus complexe que 
celle qu’on leur attribue généralement. Ces couches ne sont pas continues ; leur structure granulaire est mise en évidence par le 
fait qu’elles diffusent une partie de la lumiére qu’elles recoivent, ce qui peut donner lieu a un tres brillant systéme de franges de 
diffusion dans la zone de réflexion totale ; elles ne sont pas homogeénes, c’est-a-dire présentent des gradients d’indice et des couches 
de passage ; enfin, elles sont quelquefois anisotropes. ; 


ZUSAMMENFASSUNG. — Die optischen Eigenschaften aufgedampfler ditnner durchsichtiger Schichten werden experimentell unter- 


sucht. Die Ergebnisse erlauben die Struktur solcher Schichten zu bestimmen. Dabei werden spektrophotometrische Methoden bei 
senkrechtem und schiefem Lichteinfall benutzt. Besonders eingehend untersucht sind die Erscheinungen nahe der Grenze der 
Totalreflexion. 

Alle Ergebnisse deuten darauf hin, dass die untersuchten dtinnen Schichten eine sehr viel verwickeltere Struktur haben als 
man gemeinhin glaubt. Sie sind nicht homogen ; ihre kirnige Struktur geht daraus hervor, dass sie einen Teil des einfallenden 
Lichtes streuen. An der Grenze der Totalreflexion kénnen sie zu sehr hellen Diffusionsstreifen Anlass geben. Sie zeigen neben 


ungestorten Durchgdngen ortliche Verdnderungen des Brechwertes, also Anisotropie. 


SUMMARY. — An experimental study has been made of some optical properties of thin transparent films formed by evaporation, and 
it is shown that these results can be used to determine the structure of the films. Spectropholometric measurements have been made 
at normal and oblique incidence, with particular attention to the region of total reflection. 

The results show that the thin films studied have a much more complex structure than has generally been supposed. These films 
are not continuous, their granular structure giving rise to diffusion of some of the incident light, which can give rise to a brilliant 
system of diffusion fringes in the region of total reflection ; they are inhomogenous, showing refractive index gradients and “‘ cou- 


ches de passage ”’ ; finally they are sometimes anisotropic. 


Introduction. — On admet le plus souvent qu’une 
couche mince transparente déposée par évaporation 
sous vide est assimilable du point de vue optique a 
une simple lame a faces planes et paralléles, comple- 
tement définie par son indice et son épaisseur. C’est 
évidemment l’hypothése la plus simple. Elle permet de 
prévoir facilement, par des calculs d’interférences 
classiques, les diverses propriétés optiques des couches, 
et, inversement, d’imaginer de nombreuses méthodes 
pour déterminer les deux paramétres caractéristiques 
a partir de données expérimentales. Beaucoup d’au- 
teurs ont travaillé dans cette voie. 

Malheureusement, il n’a jamais été fait, A ma con- 
naissance, de vérification expérimentale précise de 
cette hypothése fondamentale. Or, il faut bien se rendre 
compte que lorsqu’on admet que la couche est parfai- 
tement définie par son indice et son épaisseur, on sup- 
pose implicitement qu’au moins les conditions sui- 
vantes sont remplies : 


1° Ja couche est continue, puisqu’on néglige la lu- 
-miére diffusée qui pourrait étre due a une structure 


granulaire ; 
20 elle est homogéne, c’est-a-dire d’indice constant 


dans toute l’épaisseur ; 
3° chaque surface de séparation, cdté air et cote 
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support, correspond a une discontinuité parfaite de 
Vindice, sans existence de couches de passage ou zones 
de transition d’aucune sorte ; 


49 la couche est tsotrope. 


Or, rien ne prouve a priori que les couches minces 
obtenues par évaporation thermique possédent toutes 
ces propriétés. Etant donné le mécanisme de leur 
préparation, cela semblerait méme assez étonnant. 
Nous nous sommes par suite proposés les buts suivants 


fpedsliges 


1° déterminer d’abord, avec la plus grande précision 
possible, les propriétés optiques d’un certain nombre 
de corps transparents en couches minces ; 


2° comparer les résultats obtenus avec ceux qu’on 
peut calculer en appliquant la théorie simplifiée de la 
lame continue, homogene et isotrope ; 


3° montrer enfin comment les différences observées‘ 
parfois trés importantes, peuvent s’expliquer si on 
abandonne les hypotheses simplificatrices habituelles, 
en déduire des indications sur la structure véritable 
des couches et prouver, a l’occasion, comment les métho- 
des classiques de détermination de l’indice et de l’épais- 
seur des lames peuvent conduire a des résultats erronés. 
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Nous diviserons cet exposé en deux parties princi- 
pales. La premiére sera consacrée aux résultats fournis 
par les mesures spectrophotométriques effectuées en 
incidence normale ou en incidence oblique. Dans la 
deuxiéme partie nous mettrons en évidence les nom- 
breux résultats que fournit l’étude des couches minces 
transparentes dans la zone de réflexion totale. 


I. Mesures spectrophotométriques. — L’appareil 
avec lequel ces mesures ont été effectuées [1] a été 
concu pour répondre a un certain nombre d’exigences : 
il permet les mesures de facteurs de réflexion et de 
transmission sous toutes les incidences depuis linci- 
dence normale jusqu’au voisinage de l’ incidence rasante ; 
le faisceau lumineux monochromatique est polarisé 
rectilignement avant de tomber sur la couche mince ; 
enfin, lincertitude sur les mesures de facteurs de 
réflexion et de transmission ne dépasse pas deux a 
trois milliemes en valeur relative, dans le cas général. 


1. Mesures en incidence normale. — Dans le cas 
idéal de la lame continue, homogene, a faces planes et 
paralléles, d’indice n, et d’épaisseur d, bien définis, 
le calcul montre facilement que les facteurs de réflexion 
R et de transmission 7’, doivent subir, en fonction de la 
longueur d’onde de la lumiére incidente, des varia- 
tions d’allure sinusoidale. C’est la un résultat classique, 
sur lequel est basée une méthode, classique également, 
de détermination des quantités n, et d,. Si Vindice n, 
de la couche a une valeur comprise entre les indices 
n, et n; des milieux environnants, on montre que les 
maxima des courbes donnant le facteur de réflexion R 
en fonction de la longueur d’onde ont pour expres- 
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cest-a-dire sont égaux au facteur de réflexion du sup- 
port nu. 

Nous avons mesuré, en fonction de la longueur 
d’onde, les facteurs de réflexion sous incidence nor- 
male d’un grand nombre de couches minces transpa- 
rentes, d’épaisseurs variées, obtenues par évaporation 
thermique de diverses substances, telles que le fluo- 
rure de calcium, la cryolithe, le fluorure de magnésium. 
Les résultats que nous avons obtenus ne s’accordent 
jamais avec ceux que laissent prévoir les calculs basés 
sur ’hypothése de la lame continue et homogéne. En 
particulier les maxima du facteur de réflexion de la 
couche ne sont jamais égaux au facteur de réflexion du 
support. 

Les exemples les plus typiques sont fournis par les 
couches de fluorure de calcium pour lesquelles les 
maxima de # sont trés inférieurs au facteur de réflexion 
du support. De plus, ces maxima ont une valeur qui 
varie avec la longueur d’onde (fig. 1 et 2). 

Sur la figure 2 on a reporté également la courbe 
représentant, en fonction de la longueur d’onde, le 
facteur de transmission 7’ de la couche. On peut véri- 
fier, qu’a la précision de nos mesures photométriques 
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(2 a3 milliémes en valeur relative), R et 7 sont par- 
faitement complémentaires ; il n’y a done pas d’absorp- 
tion ni de diffusion mesurables pour de telles couches, 
en incidence normale. Ceci a une conséquence pratique 


4 
R 
go44 
90034 
goe2 | 
Oo14 
AA 
Ou SR te © 5 T => 
750) 6150 5100 4485 4055 3720 
Fic. 1. — Facteurs de réflexion en fonction de la 


longueur d’onde d’une couche de fluorure de calcium, 
d’épaisseur optique voisine de 9 600 A, déposée sur 
verre. 


importante : par suite de la faible valeur des maxima 
de R et de absence d’absorption, les minima de 7 sont 
plus grands, de 2/100 environ, que ne le laisse croire la 
théorie classique. L’effet « antireflet » de ces couches de 
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Fic. 2, — Facteurs de réflexion et facteurs de trans- 
mission, en fonction de la longueur d’onde, d’une 
couche de fluorure de calcium, d’épaisseur optique 
voisine de 6 900 A, déposée sur quartz. En pointillé : 
facteurs de réflexion Rg et de transmission Tq du 


support. 


fluorure de calcium est donc bien meilleur qu’on ne le 
crout habituellement. La figure 3 représente, en fonction 
de la longueur d’onde, les variations du facteur de 
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réflexion d’une couche de cryolithe et la figure 4 celles 
du facteur de réflexion d’une couche de fluorure de 
magnésium. On constatera que les maxima de R sont, 
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Fic. 3. — Rapport du facteur de réflexion de la 
couche mince a celui du support, en fonction de la 
longueur d’onde, pour une couche de cryolithe,d’épais- 
seur optique voisine de 8 000 A, déposée sur quartz. 


la encore, différents du facteur de réflexion du support 
nu ; mais, contrairement au cas du fluorure de calcium, 
ils lui sont ici supérieurs. Il convient d’ailleurs de 
noter que, pour le fluorure de magnésium, l’écart entre 
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Fic. 4. — Rapport du facteur de réflexion de la 
couche mince a celui du support, en fonction de la 
longueur d’onde, pour une couche de fluorure de 
magnésium déposée sur veire. 


es maxima du facteur de réflexion de la couche et celui 
du support est beaucoup plus faible que pour la cryo- 
lithe et le fluorure de calcium. 

Les résultats expérimentaux que nous venons d’expo- 
ser sont inexplicables dans lhypothése d’une couche 


PROPRIETES OPTIQUES DES COUCHES MINCES TRANSPARENTES 155 


homogene d’indice n, constant. Ils peuvent par contre 
s’interpréter correctement si l’on abandonne cette 
hypothése de Phomogénéité des couches minces. 

L’existence d’un simple gradient d’indice, c’est-a- 
dire d’une variation réguliére de n, d’une face a l’autre 
dela couche, pourrait suffire 4 expliquer le fait que les 
maxima du facteur de réflexion sont différents du 
facteur de réflexion du support. Mais cette hypothése 
ne peut expliquer la variation de la grandeur de ces 
maxima en fonction de la longueur d’onde ; cette 
variation est pourtant trés importante. L’ensemble 
des résultats expérimentaux rappelés ci-dessus peut 
par contre étre expliqué completement si lon admet 
Pexistence de couches de passage sur les surfaces de 
séparation de la couche mince avec les deux milieux 
qui Pentourent. Nous verrons d’ailleurs plus loin que 
les résultats des mesures spectrophotométriques en 
incidence oblique confirment, de facon plus directe et 
sans aucun doute possible, existence de ces couches 
de passage. On peut alors montrer que si 9’ el? et o” el? 
désignent les coefficients de réflexion complexes de 
chacune des couches de passage, pris du cété de la 
couche principale, le facteur de reflexion de ensemble 
s’éerit : 
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Si les couches de passage n’ont pas une épaisseur tres 
petite par rapport a la longueur d’onde, .’ et pe” 
varient avec A et il s’en suit que les maxima de RF sont 
également fonction de 4. 
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Il est possible, tout au moins si l’épaisseur de la 
couche est suffisante, de déduire de la détermination 
expérimentale des facteurs de réflexion des renseigne- 
ments suffisamment précis sur les caractéristiques de 
la couche principale et des couches de passage. La 
figure 5 résume le résultat de ces déterminations pour 
la couche de fluorure de calcium dont les facteurs de 
réflexion sont représentés par la courbe de la figure 1. 

On a la un exemple des résultats que l’on peut tirer 
de l’étude spectrophotométrique des couches minces 
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transparentes, en incidence normale. II faut cependant 
remarquer que ces résultats ne sont pas complets, car 
nous avons supposé, pour faire la figure ci-dessus, que 
Vindice de la couche principale était constant ; or, 
rien ne prouve qu’il en est ainsi et il est fort possible 
que, comme la couche de passage, cette couche prin- 
cipale présente elle-méme un gradient d’indice. 

Remarquons, d’autre part, que nous sommes ici 
dans un cas relativement simple, puisque nous avons 
une seule couche de passage, sur Vinterface air-fluo- 
rure. Sur l’interface fluorure-verre, la couche de pas- 
sage, si elle existe, est d’épaisseur tres petite et par 
suite, n’a pas d’influence sur la valeur des facteurs de 
réflexion. I] n’en est pas toujours ainsi et il arrive que 
les deux couches de passage existent et aient toutes 
deux une influence non négligeable. 

A la lumiére des résultats fournis par les mesures de 
facteurs de réflexion sous incidence normale, on voit 
que les couches se présentent généralement comme un 
ensemble complexe constitué par une couche prin- 
cipale, dont l’indice n’est pas obligatoirement constant, 
entouré d’une ou de deux couches de passage, elles- 
mémes non homogenes. 

Méme dans le cas favorable de couches assez épaisses 
il est impossible de déduire du résultat de ces seules 
mesures tous les renseignements nécessaires sur la 
structure de ces couches. C’est que, par suite de la 
complexité de cette structure, le nombre de paramétres 
a déterminer devient important. Nous sommes loin 
de la simple lame homogéne, parfaitement définie par 
deux grandeurs seulement, l’indice et Vépaisseur. I] 
s’en suit que l’application brutale a de telles couches 
des méthodes classiques de détermination des para- 
métres caractéristiques ne peut conduire qu’a des 
résultats erronés. Nous allons voir maintenant com- 
ment les résultats des mesures spectrophotométriques 
en incidence oblique viennent confirmer et compléter 
ceux fournis par les mesures en incidence normale. 


2. Mesures en incidence oblique. — Comme en inci- 
dence normale, il est facile de montrer que pour une 
incidence 6, quelconque, le facteur de réflexion R 
d’une couche mince transparente, homogéne et iso- 
trope, subit, en fonction de la longueur d’onde, des 
variations d’allure sinusoidale. Supposons encore l’in- 
dice n, de la couche intermédiaire entre les indices 
des milieux environnants. Sila vibration lumineuse est 
perpendiculaire au plan d’incidence, les maxima de R 
sont égaux au facteur de réflexion du support. I] en 
est de méme dans le cas ot la vibration est paralléle 
au plan d’incidence, si l’angle d’incidence 0, est infé- 
rieur a langle de BREwsTER 0, sur la surface air-couche 
mince. Pour les incidences supérieures, ce sont les 
minima de A) qui sont égaux au facteur de réflexion 
du support. 

Les mesures spectrophotométriques en incidence 
oblique pourraient permettre l’étude de l’anisotropie 
dune couche mince homogéne. Nous avons en effet 
montré | 1] qu’a partir du résultat de ces mesures de 
facteurs de réflexion, on peut calculer, pour n’importe 
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quelle incidence, les indices de la couche relatifs a ses 
différents axes de symétrie électrique. 

Malheureusement, une telle méthode n’est appli- 
cable que dans le cas d’une couche homogéne. Nous 
savons déja que, pratiquement, il n’en est généralement 
pas ainsi ; ici aussi l’interprétation des résultats n’est 
donc pas aussi immédiate que ne le laissent croire les 
formules. 

Nous avons effectué des mesures de facteurs de 
réflexion en incidence oblique principalement sur des 
couches de cryolithe et de fluorure de calcium. De ces 
mesures nous pouvons tirer trois résultats principaux 
qui sont en contradiction avec ’hypothése dela couche 
homogene et isotrope. 


a) Comme en incidence normale, les maxima du 
facteur de réflexion R de la couche (ou les minima s’il 
s’agit de la vibration paralléle et d’incidences supé- 
rieures a 6,), ne sont jamais égaux a leur valeur théo- 
rique qui est le facteur de réflexion du support. 
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Fic. 6. — Courbes théoriques donnant, en fonction 
de la longueur d’onde, le rapport R r/Re du facteur 


de réflexion d’une couche mince a celui du support 
nu, pour diverses incidences voisines de Vincidence 
brewstérienne sur la surface air-couche mince. Ces 
courbes ont été calculées dans le cas d’une couche 
homogéne Windice 1,375. 


b) Il n’y a pas disparition des interférences dans la 
lame pour la vibration paralléle au plan d’incidence 
lorsque l’angle d’incidence est égalal’angle de Brews- 
TER sur la surface air-lame [14 ] [ 2]. 

c) Les positions des maxima et des minima de R 
dans le spectre sont, en général, différentes pour les 
deux directions de Nolan e tion [a ee 


Le premier de ces résultats est normal ; il n’est que 
Pextension de ce que nous avons déja vu en incidence 
normale. 

Le deuxiéme résultat apparait comme une preuve 
directe de lexistence d’une couche de passage sur la 
surface de séparation lame-air. I] est facile de montrer 
en effet que, dans le cas d’une couche homogéne, les 
interférences dans la lame doivent disparaitre lorsque 
0, est égal a 6,, le coefficient de réflexion sur l’interface 
lame-air étant nul pour la vibration paralléle. Dans ces 
conditions, le facteur de réflexion de la couche mince 
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Rj, sous Pincidence 6,, doit étre égal au facteur de ré- 
flexion du support pour toutes les longueurs d’onde [5}. 
Les courbes de la figure 6 représentent les variations 
du rapport R,/R,. du facteur de réflexion de la couche 
mince a celui du support, calculées pour quelques 
angles d’incidence au voisinage de 0,, dans le cas d’une 


couche mince supposée homogene, d’indice 1,375 
(8, = 54°). On a tenu compte, dans le calcul, de la 
dispersion de la couche mince, ce qui explique les petits 
écarts entre la courbe relative a 6, = 54° et la droite 
d@ordonnée 1. On voit que, dans le cas d’une lame par- 


R 
fartement homogene, les oscillations du rapport oT en 


Q 
fonction de la longueur d’onde doivent s’effacer sur 


place 4 mesure que 0, croit. jusqu’a 6,, pour s’inverser 
lorsque 0, est supérieur a 6,. 
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Fic. 7.— Rapport du facteur de réflexion d’une cou- 
che de cryolithe 4 celui du support nu, en fonction 
de la longueur d’onde, pour différentes incidences 
voisines de l’incidence brewstérienne sur la surface 
air-couche mince (courbes expérimentales). 


Voyons ce qui se passe en réalité : sur la figure 7 on 
a tracé les courbes expérimentales représentant le 
rapport du facteur de réflexion de la couche mince a 
celui du support, en fonction de Ja longueur d’onde, 
pour différentes incidences voisines de lincidence 
brewstérienne sur la surface air-couche mince. Le 
rapprochement de ces courbes expérimentales avec les 
courbes de la figure 6 montre que le passage par l’inci- 
dence 0, se fait de facon tout a fait différente de celle 
attendue pour une couche homogene. Au lieu d’avoir 
passage d’un maximum a un minimum, ou inverse- 
ment, par déformation sur place de la courbe, on a en 
réalité déplacement de l’ensemble de cette courbe vers 
les courtes longueurs d’onde. Les oscillations du rap- 
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J . . 
port Rene disparaissent pour aucune valeur de l’angle 


Q 
d@incidence, mais restent, au contraire, toujours tres 
marquees. 

Ce résultat est sans doute la preuve la plus directe 
de l’existence d’une couche de passage sur l’interface 
lame-lair. Le fait que les oscillations du facteur de 
réflexion Ry ne disparaissent jamais prouve que le 
coefficient de Fresnev sur linterface air-lame n’est 
nul pour aucune incidence ; le déplacement progres- 
sif des maxima et des minima de R, dans le spectre 
prouve que le changement de phase a la réflexion sur 
cet interface varie progressivement au lieu de présen- 
ter une discontinuité de =. Ce sont la des résultats 
caractéristiques de l’existence d’une couche de pas- 
sage sur la surface de séparation considérée. 
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Fic. 8. — Positions dans le spectre, en fonction de 
angle d’incidence et pour les deux directions prin- 
cipales de polarisaticn, des maxima et minima du fac- 
teur de réflexion d’une couche mince de fluorure 
de calcium. 


Cette existence est confirmée par notre troisieme 
résultat expérimental : les positions des maxima et 
minima du facteur de réflexion dans le spectre sont 
différentes pour les deux directions de polarisation. Ceci 
est illustré par les courbes des figures 8 et 9. La pre- 
miére représente les positions dans le spectre, en fonc- 
tion de l’angle d’incidence, des maxima et des minima 
du facteur de réflexion d’une couche mince de fluo- 
rure de calcium, pour les deux directions principales 
de polarisation. La figure 9 donne les mémes résultats 
pour une couche de cryolithe. On retrouve sur ces 
courbes le phénomene déja décrit du déplacement 
rapide dans le spectre des maxima et minima de Rj 
au passage par lincidence 0,. On y voit, de plus, qwil 
existe en général pour toutes les incidences un écart 
Ad plus ou moins grand entre les positions des maxima 
et des minima relatifs aux deux directions principales 
de polarisation. Cet écart Ad est de signe différent 
suivant que l’angle d’incidence est inférieur ou supé- 
rieur a 0,. Ce dernier résultat exclut toute possibilité 
dinterpréter ce résultat uniquement en faisant appel 
a une anisotropie possible des couches. Par contre, 
Yensemble des résultats représentés par les courbes 
ci-dessus s’explique complétement a partir de l’exis- 
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tence d'une couche de passage sur linterface couche- 
air. Les écarts AX observés entre les positions des 
maxima et des minima relatifs aux deux directions de 
polarisation apparaissent alors dus aux différences 
entre les changements de phase a la réflexion sur la 
couche de passage relatifs aux deux composantes. 
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Fic. 9. — Positions dans le spectre, en fonction 
de langle d’incidence, et pour les deux directions 
principales de polarisation, des maxima et minima 
du facteur de réflexion d’une couche mince de cryo- 
lithe. 


Les différents résultats que Von peut tirer des 
mesures spectrophotométriques en incidence oblique 
confirment done ceux que l’on avait déduit des mesures 
en incidence normale. Nous avons pu vérifier accord 
entre les valeurs numériques des paramétres carac- 
téristiques des couches étudiées,qu’elles soient déduites 
du résultat de mesures en incidence normale ou en 
incidence oblique. Ces deux types de mesures peuvent 
se compléter mutuellement et nous permettre d’ac- 
croitre nos connaissances sur la structure des couches 
minces transparentes. 

Remarquons enfin que les mesures spectrophoto- 
métriques ou polarimétriques en incidence oblique 
sont ala base de plusieurs méthodes de détermination 
de l’indice et de l’épaisseur des lames minces transpa- 
rentes. Toutes ces méthodes sont naturellement basées 
sur Vhypothése que la couche étudiée est homogeéne 
et isotrope. A la lumiére de ce que nous savons main- 
tenant de la structure de ces couches, ces méthodes 
apparaissent d’application difficile ; elles peuvent 
souvent conduire a des erreurs. 


Il. Etude des couches minces transparentes dans la 
zone de réflexion totale. — Les couches minces n’ont 
été que trés rarement étudiées sous des incidences 
sipérieures a l’angle limite, sauf peut-étre en ce qui 
concerne le cas de la réflexion totale « frustrée » qui 
a trouvé des applications dans les filtres interférentiels. 
Nous allons montrer que l’on peut pourtant tirer des 
résultats intéressants d’études faites dans cette région. 


BOUSQUET 
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1. Premiere Experience Franges par défaut de 
réflerion totale. —- Considérons Vexpérience schéma- 
Lisée par la figure 10. Une couche de fluorure de cal- 
cium, d’épaisseur supérieure a 2 % environ, est déposee 
sur la face hypoténuse d’un prisme rectangle isocele. 
Un faisceau lumineux convergent, polarisé rectiligne- 
ment, parallélement ou perpendiculairement au plan 
d’incidence, se réfléchit sur cette face. I est facile de 
montrer que, théoriquement, la présence d’une couche 
mince transparente sur Phypoténuse du prisme ne 
doit en rien modifier le phénomene observé sur |’écran. 
La séparatrice délimitant la zone de réflexion totale 
de celle de réflexion partielle doit occuper la méme 
place que dans le cas du prisme nu ; au dela de cette 
séparatrice la réflexion doit étre totale. 


Poloriseur S, 


Ecran 
Fie. 10: 


En réalité, si Vhypoténuse du prisme est recouverte 
dune couche de fluorure de calcium d’épaisseur suffi- 
sante, la réflexion n’est pas totale au dela de la sépa- 
ratrice. Tout une région de l’écran apparait sillonnée 
de franges sombres bien visibles. On peut vérifier que 
cest la région correspondant aux angles d’incidence 
compris entre Vincidence limite fluorure-air et Vinci- 
dence limite verre-fluorure. Au dela de cette derniere 
incidence limite, la réflexion est effectivement totale. 

La figure 11 est la photographie du systeme de 
franges obtenu avec une couche de fluorure de calcium 
d’environ 10 A d’épaisseur optique (4 = 4358 A, vibra- 
tion perpendiculaire au plan d’incidence). Des mesures 
de facteurs de réflexion que nous avons effectuées 
dans cette zone précisent quantitativement ce phé- 
nomene. Au lieu d’étre égal a 1 pour toutes les inci- 
dences supérieures 4 l’incidence limite couche-air, le 
facteur de réflexion présente des minima trés nets, 
d’autant plus nombreux et mieux marqués que I’épais- 
seur de la couche est plus grande. 

Ce phénomene trouve son explication dans le fait 
que, par suite de leur structure granulaire, les couches 
minces, particuliérement celles de fluorure de calcium, 
diffusent une partie de la lumiére qui les traverse. C’est 
a cette perte de lumiére par diffusion qu’est dt le fait 
que le facteur de réflexion, au dela de langle limite, 
nest pas égal 4 1. Nous avons pu montrer que l’on 
retrouve par le calcul Vallure des courbes expéri- 
mentales donnant le facteur de réflexion au dela de 
Yangle limite si ’on affecte a la couche un indice 
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Fig. 11. — Franges sombres obtenues par réflexion sur une couche de fluorure de calcium, 
pour des incidences supérieures a angle limite (expérience de la figure 10). 


d’extinction » pour traduire la perte de lumiére due a 
la diffusion. Cet indice d’extinction est toujours tres 
faible (de ordre de 10-4 environ), ce qui explique quwil 
ne produise aucun effet mesurable en incidence nor- 
male. Par contre, son importance devient prépondé- 
rante ici, du fait de l’inclinaison beaucoup plus grande 
des rayons dans la lame. C’est donc la exemple d’un 
phénomene que seules des mesures dans la zone de 
réflexion totale permettent d’étudier facilement. 

La lumiere diffusée par les couches minces de fluo- 
rure de calcium peut étre mise en évidence de fagon 
particuliérement brillante par une autre expérience, 


Fig. 12. 


schématisée par la figure 12 [4]. La couche mince, 
toujours déposée sur la face hypoténuse du prisme, 
est éclairée par un faisceau lumineux d’ouverture 
et d’incidence quelconques. On s’arrange simplement 
pour qu’il ne puisse pas y avoir de rayon qui subisse 
la réflexion totale sur la face hypoténuse. Dans ces 
conditions, en l’absence de lumiere diffusée, une lunette 
visant dans des directions situées dans la zone de 
réflexion totale, ne doit recevoir aucune lumiére. En 
réalité, au lieu d’étre obscur, le champ dans la lunette 
apparait sillonné de franges brillantes, fines et tres 


contrastées. On peut vérifier que ces franges brillantes 
de diffusion occupent exacternent la place des franges 
sombres observées dans la lumiére réfléchie lors de la 
premiere expérience. Ces franges de diffusion sont dues 
au fait que la lumiere diffusée par chaque particule de 
matiére subit des réflexions multiples dans la couche et 
se répartit par suite suivant un systeme de franges a 
Pinfini. Observées en lumiere blanche, ces franges pré- 
sentent de vives colorations. La figure 13 représente 
les franges de diffusion produites par la couche de 
fluorure de calcium ayant fourni les franges par défaut, 
de réflexion totale de la figure 11; la longueur d’onde 
et état de polarisation de la lumiére sont également 
les mémes. On peut vérifier facilement que les deux 
systemes sont complémentaires. 


2. Deuxieme Experience. Anomalies au voisinage de 
Vangle limite verre-couche : Franges par polarisation. — 
Les phénomenes que nous venons de décrire s’inter- 
prétent correctement, en tenant compte de la diffusion 
de la lumieére a Vintérieur de la couche, sauf pour les 
incidences les plus grandes, c’est-a-dire au voisinage 
de l’angle limite verre-fluorure. Dans une zone d’un ou 
deux degrés avant cette incidence limite, on observe 
des irrégularités dans la position des franges, dont la 
distribution ne peut plus s’expliquer aussi simplement 
que nous l’avons fait jusquw’ici. Une nouvelle expérience 
va nous prouver l’existence de phénomenes plus 
complexes dans cette région. 

Nous reprenons le montage de notre premiere expé- 
rience (fig. 10), mais nous plagons maintenant deux 
polariseurs de part et d’autre du prisme, leur section 
principale étant a 45° du plan d’incidence. Le calcul 
laisse prévoir que nous obtiendrons sur l’écran un 
systeme de franges di a des interférences en lumiére 
polarisée. Mais ces franges ne devraient avoir qu’un 
contraste faible ; leur répartition dans le plan d’obser- 
vation devrait d’autre part étre réguliere. Or, on observe 
en réalité un systeme de franges excessivement con- 
trasté présentant des minima nuls et des maxima tres 
intenses. En lumiére blanche le phénoméne est parti- 
culiérement brillant. De plus, au voisinage de langle 
limite verre-fluorure, la position de ces franges n’est 


Fic. 13. — Franges de diffusion obtenues avec la méme couche de fluorure de calcium 
que dans le casde la figure 11 (expérience de la figure 12). 
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pas du tout celle prévue par la théorie simple faite 
jusqu ici. 

La répartition de l’intensité lumineuse dans le plan 
d’observation au cours de lexpérience précédente est 
conditionnée par la valeur de la différence Ar, — Ar) 
entre les changements de phase par réflexion sur la 
couche mince relatifs aux deux directions principales 
de polarisation. Il] semble done que, dans la région 
voisine de l’incidence limite verre-fluorure, cette diffé- 
rence Ar, — Ar) subisse des variations anormales 
et beaucoup plus importantes que celles prévues par le 
calcul. Ceci a été confirmé par des mesures de la quan- 
tité Ar, — Ar réalisées a Taide d’un analyseur de 
CuAumonr. Les mesures ont montré effectivement que 
lorsque Vépaisseur de la couche dépasse 1 ou 2 4, il 
y a toujours au voisinage du 2° angle limite, un désac- 
cord important entre les valeurs de Ar, — Ar mesu- 
rées et celles que l’on peut calculer a partir des données 
que l'on posséde sur la couche ¢tudiée. Les courbes 
experimentales présentent, en particulier, d’impor- 
tantes oscillations que l’on ne retrouve jamais sur les 
courbes théoriques. 

Rien de ce que nous avons appris jJusqu’ici de la 
structure de nos couches ne nous permet d’expliquer 
ces anomalies. Il semble que seule ’hypothese d’une 
anisotropie de la couche mince puisse fournir l’expli- 
cation cherchée. Nous avons pu montrer en effet qu il 
suffisait de supposer l’existence d’une tres faible ani- 
sotropie de ja couche (différence d’indices de 3 a 4 mil- 
liemes environ) pour retrouver, par le calcul, exacte- 
ment la forme des courbes expérimentales représen- 
tant Ar, —Ar) en fonction de Pangle d’incidence. 


Il en est done de l’anisotropie comme de la diffu- 


[ Orr. AcTA 
sion ; inappréciables lune et l'autre en incidence nor 
male ou en incidence oblique ordinaire, elles prennent 
en réflexion totale une importance prépondérante. 

Par suite de la trés grande inclinaison des rayons 
lumineux dans la couche mince toute une série de 
phénomeénes trop peu importants pour étre décelables 
dans les conditions ordinaires prennent ici une impor- 
tance considérable. C’est dire tout Vintérét qui s’at- 
tache a des mesures effectuées dans cette zone. 


Conclusion. — Les résultats de toutes les mesures 
que l’on peut effectuer sur les couches minces trans- 
parentes concourent a prouver que la structure de ces 
couches est bien plus complexe qu’on ne ’admet géné- 
ralement. Elles ne sont pas homogenes, c’est-a-dire 
présentent des gradients d’indice et des couches de 
passage ; elles ne sont pas continues, leur structure 
granulaire ¢tant mise en évidence par le fait qu’elles 
diffusent une partie de la lumiere qu’elles recoivent ; 
enfin, elles apparaissent quelquefois anisotropes. Par 
suite, le nombre de parametres qui les caractérise est 
tres grand ; les méthodes classiques de détermination 
des constantes optiques de ces couches sont alors le 
plus souvent en défaut. 


REFERENCES 


1] P. Bousquet, Thése Doct. Sc. Phys. Paris, 1956 (a parai- 
tre aux Annales de Physique). 

] P. Bousquet, C. R. A. S., 240, 1955, p. 2502. 

] P. Bousquet, C. R. A. S., 241, 1955, p. 478. 

] P. Bousquet, C. R. A. S., 287, 1953, p. 516. 

] F. ABEvEs, C. R. A. S., 228, 1949, p. 553. 


Manuscrit recu le 16 juillet 1956. 


vol. 8, n° 4 Dec. 1956] TOTAL ILLUMINATION IN AN ABERRATION-FREE DIFFRACTION IMAGE 161 
Total illumination in an aberration-free diffraction image 
J. Focke 
Mathematisches Institut, Leipzig 
SUMMARY. — Asymptotic expressions are derived for the total illumination in a defocused aberration-free diffraction image based 
on the asymptotic approximations of the diffraction integral due to SCHWARZSCHILD or 10 VAN KAMPEN respectively. 
SoMMAIRE. — Des expressions asymptotiques basées sur les approximations de Vintégrale de diffraction de ScHWARZSCHILD ou de 


VAN KAMPEN sont données ici pour représenter énergie « encerclée » dans une image de diffraction, dépourvue d’aberrations, 


en présence d’un défaut de mise au point. 


ZUSAMMENFASSUNG. — Es werden auf Grund der von SCHWARZSCHILD bzw. YAN KAMPEN angegebenen asymptotischen Ndherungen 
des Beugungsintegrals asymptotische Ausdritcke fiir die Totalbeleuchtung in einem extrafokalen aberrationsfreien Beugungsbild 


hergeleitet. 


Beside the intensity distribution in an aberration- 
free diffraction image, which was studied in classical 
memoirs by Airy (1835) and Lommet (1885) and more 
recently by ZernIkE and NiJBOER, it is also of interest 
to know the total illumination in the various rings of 
the diffraction pattern. In modern times formulae 
were published by E. Wotr [1] which permit direct 
calculation of the total illumination. In the present 
paper approximations are derived for the total illumi- 
nation, which represent the first term of asymptotic 
expansions with respect to small wavelength. These 
approximations therefore are intermediate between 
physical and geometrical optics. They asymptotically 
approach the exact results of wave theory and can 
often be used instead of them. The numerical compu- 
tations are of a very simple kind, even at large dis- 
tances from the focus. 

Following [1] we introduce the following notation : 


f focal length, 

Af distance between the receiving plane and the 
geometrical focus, 

R_ radius of the opening of the aperture, 

r, 8 polar co-ordinates in the receiving plane, 

2 wave-length, 

u, 9 normalized co-ordinates, 


Dh ae 18: 
“21, (0 =0) 


iN 


(u = v, boundary of the geometrical shadow), 


4 (2 (2R . z 
<s | | e—ilve cos p+up2/2) 0 do de 
Te) (01 9/0 


light intensity, normalized to )(0,0) = 1 ; 


(1) Lu, ¢) = xf. (u,v) 9 de 


total illumination of a disk with radius Vo 
in the receiving plane wu = const., normalized 
to L(u, ©)=1. 


J(u, ¢) = 


We start from the asymptotic expansions of the 
diffraction integral with respect to wu and ¢ (impli- 
citly containing 1/2) which were first published by 
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N. G. van Kampen [ 2] and later on by the author [3] 
im a rigorous mathematical treatment. Under re- 
striction to the first term we have 


4 
2) °3@,%)~—, for sabe 

uU 

9 lie ile—7/4) j ell—7/4) |? 
(3) ju, 9) y — 

Te g+u 9—wU 


~e| | 4 4 : 4sin2¢ 
my} (e+u)? (vy —u)?  -79(92—u?) 


for 0<u<¢. 


These expressions in the error-free case are also given 
by ScHwarzscHILp [4]. They fail to be valid in the 
neighbourhood of the edge of the geometrical shadow. 
According to [2] or to [4] in that region the intensity 
can be approximated by FrESNEv’s integrals. 


kel 2 teria! 


g+u 


, ew? = 
where B=e 2 4 x 
dl ‘ 5 fal 
NK | ( 5 CH) ) i . 


t 
Ci) + i S(t) = | ent qe, 


2 


(4) Su, ¢) + Beil*—n/) 


ne | 
forr0<u<e 


This formula not only is valid near the edge of the 
geometrical shadow but also in the whole shaded region, 
for by substituting 


(5) IG = CA)! (5 = so) | ee 


7 : 
i , 
— — eine fort >0 
~ Tt 


1—i fort <0 


expression (4) passes into (3) when ¢ > u. Within the 
illuminated region we have similarly 


i) 
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(6) 27/4 \2 4 : 
wu) ~ 2 (= —c()) + (5-50) A emu 


Vi Tu 


for0<e<u 


and this expression is also valid at a large distance from 
the edge of the shadow ; for applying (5) formula (4) 
goes over into (2) when »<u. On the boundary of the 
shadow the first term from (4) is of a smaller order of 
magnitude than the second term and can be neglec- 
ted. Equations (4) and (6) are therefore in good agree- 
ment for uv = ». However all approximations con- 
structed above fail to be true in the neighbourhood of 
the optical axis. Consequently we calculate the total 
illumination in the following form : 


(7) Llu, ©) — 


The case u < . — ae ae from (3) into (7) 
we obtain 
1 E 1 ra 1 
Tt ie ae Fil 


sh sin 2 ¢ 
By partial integration 


sine | => COS 2 V5 cos 2 9 
@ f =e = y dy. 
bo are 2 9 —u? 


We find that the integral in (8) is of a smaller order 
of magnitude than the first term and can therefore 
be neglected. We finally obtain 


(8) 1— L(y, %) ~ 


re re 


1— L(u, %) ~ ogeels for 0< 


2 772 
Tv oe u 


(10) 


SUV, & 


The case u < Vo, (Vo — u) small. — Substituting 
from (4) into integral (7) we obtain, to the first order, 


(41) 1—L(u, ¢) ~ +f. |Bede=\/2 [x Om teu 
forO<u<%, (% —w) small ae 
where 
nf [(—08)* (sas 
== 550) eo 85 =(5—C()) sin Se 


co 
Geo) +0) 
H(0) — H(t) ~ ee = 


HG) i =r 


fort > 0, 


ire <= (Oy 


Herein the expression for H(t) is constructed by a 
double partial integration. In table 1 are given some 
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numerical values of this function. On the boundary of 
the shadow we have 


1 


TU 


(12) {— L(u, u) ~ 


The case ¥) <u. — In this case we split integral (7) 
into 


(143) 1—L(u, %) =1—L(u,u) + | 


eo 0 


Substituting from (6) into (13) we obtain 


ag 1m [y/EmO—2 70], _ os 


ru 


J(u,'e) ode. 


Vg—Uu 
Vmu 
for 0< <4, 
where 
*t 4 2 { 2 
F(t) = | (> ai c(9)) ai (F na S(:)) | de 
JOS 2 *y 
Fitj)~e fort < 0. 


Some numerical values of F(t), calculated by numeri- 
cal quadrature, are given in table 1. 


TABLE 1 
t H(t) H(—t) 
0.0 0.32 0.32 
2. 24 44 
4 18 62 
6 14 88 
8 12 1.24 
1.0 0.10 1.70 
2 08 2.24 
4 07 2.17 
6 06 oo24 
8 06 3.60 
2.0 0.05 3.91 
| 


The case 0 < 9 < U, %) small. — We have not deri- 
ved the total illumination directly from (41), because 
the asymptotic approximations fail in the neighbour- 
hood of the optical axis. But now Jj is bounded near 
the axis, therefore jv tends to zero as ¢ vanishes. 
Hence we may expect that the neighbourhood of the 
axis will not contribute essentially to the value of the 
integral (1). Consequently it is permissible to substi- 
tute from (2) into (1) and we obtain, in agreement 
with geometrical optics, 

7 


(15) Lu, 9) ~ (*:) for O0< <u. 
With the help of formulae (40), (11), (14) and (45) 
we can calculate the asymptotic values of the total 
illumination in every receiving plane wu = const., 
apart from the immediate neighbourhood of the focus. 
Naturally a distinction of several cases is necessary. 
Now we will show that the various expressions indeed 
pass into one another in the intermediary regions. 
Near the edge of the shadow for decreasing (9 — uw) 
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Fic. 1. — The total illumination L (u, v9) in the receiving planes u = 0, u = 27,u = 47. 


the term 1/(%) — uw) becomes large compared with 
1/(% + wu), therefore formula (10) is approximately 
given by a %) —u). On the other hand for increa- 
sing (¢ —u) ort ah he we have 


1 
i tb Oh epee 


and the expression (10) and (11) pass into each other. 
On the boundary of ne shadow (11) and (14) are in 
exact agreement. Within the illuminated region for 
negative increasing (%) — w) (14) passes into (45). 


(16) [Zao — Smo ny/e 2) = 


For numerical computation every expression is to be 

used until it passes into the following expression. 
Now we will make two comparisons with Wo .r’s 

results. From [4 ], (2.30) it follows that exactly 


2 =u s 2s 
(17) L(w,%)=1— & S(S) Qala 


for u < %, 


qe [ Jol”) Wea Pie 
#5 Jo+3(?) Woes): 


Substituting the asymptotic expansions of BEssEL 
functions into (17) we find that 


2s +4 
Q,,(?) , (— 1)? A(2s + 41) f 


TY 


where Q,,(¢) 


(18) 


exact behaviour (after E. Woxr [1], fig. 4a). 


asymptoticb ehaviour. 


Hence 


(19) L(w,%) ~4 5 ( u ‘s 2 ae 2 (18 
gO TV T vee 
in agreement with (10). On the edge of the shadow 
Wotr’s solution reduces to the simple form [1], 
(2.27), 


1 


(20) L(u, u) = 1—},(u) cos u —J,(u) sinu~ — =A 
in asymptotic agreement with (42). 

The asymptotic behaviour of the total illumination 
in the planes u = 0, 2x and 4 x in comparison with 
the exact behaviour due to E. Wo tr is shown in 
figure 1. 

The exact values approach closely the asymptotes 
in an oscillating manner ; the asymptotes reproduce 
therefore very well the behaviour of the total illumi- 
nation and hence the distribution of light intensity. 
In many applications, for instance for the study of the 
intensity distribution in the image of an edge, the use 
of our smooth asymptotes is more convenient than 
the use of the exact values of the total illumination. 
In comparison with geometrical optics we have obtai- 
ned a much better approximation to physical reality. 
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SUMMARY. 
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The third order aberrations in the interferometric measurement 
of concentration gradients 


Harry SVENSSON 
Laboratories of LKB-Produkter Fabriksaktiebolag, Stockholm, Sweden 


The optical path difference between two coherent rays in the interferometric recording of a diffusion boundary is com- 
puted to the accuracy represented by terms of the third degree in the entrance angle and the refractive index difference across the 
boundary. It is found that even third-degree terms acquire a minimum value for that defocusing of the middle of the cell which in 
the second order theory was found to eliminate the WiENER skewness aberration. The third order aberrations are found to be 
highly dependent on the entrance angle, and there is a certain entrance angle that is more favourable than any other angle. The 
third-order aberrations are shown to decrease rapidly with time, and it is possible to define a critical initial time period of 
diffusion, after which the aberrations can no longer be observed. This critical time period is shown to be at a minimum for the 
optimum entrance angle. Exact equations for the optimum entrance angle and for the critical time period are given in terms of 
apparatus constants and the total number of fringes in the interferogram. The blurring of fringes due to a vertically extended 
light source in conjunction with a defocused cell is put in relation to the critical time period. It is possible to arrange an auto- 
matic guard against the evaluation of interferograms not yet free from aberrations by certain adjustments causing the blurring 
and the aberration to fade out to insignificance at the same time. Finally, the effect of uncollimated light is analysed, and exact 
conditions for the divergence or convergence of the light are presented. 


Sommaire. — La différence de chemin optique entre deux rayons cohérents dans lVenregistrement interférométrique d’une frontiére 


de diffusion est calculée avec la précision représentée par les termes du troisiéme degré par rapport a Vangle d’entrée et a la diffé- 
rence d’indice de réfraction au passage de la frontiére. On trouve que les termes du troisiéme degré acquiérent une valeur minimum 
lorsque le défaut de mise au point au milieu de Vobjet a la valeur trouvée dans la théorie du second ordre comme éliminant l’aber- 
ration d’obliquité de Wiener. On trouve que les aberrations du troisiéme ordre dépendent fortement de Vangle d’entrée, et qu'il 
existe un certain angle d’entrée qui est plus favorable que tout autre. On montre que l’interférogramme est pratiquement exempt 
d aberration apres le passage d'une certaine période de temps initiale. Si Vangle d’entrée le plus favorable est utilisé, cette 
période de temps est réduite au minimum. Des équations exactes pour langle d’entrée optimum et pour la période de temps cri- 
tique sont données en fonction des constantes de Vappareil et du nombre total de franges de Vinterférogramme. Le brouillage 
des franges di a une source lumineuse étendue dans le sens vertical et a un défaut de mise au point est mis en relation avec la 
période de temps critique. Il est possible d’arranger une protection automatique contre l'emploi d’interférogrammes qui ne sont 
pas encore exempts d’aberrations, par certains dispositifs ayant pour effet que le brouillage et l’aberration deviennent négli- 
geables en méme temps. Finalement, on analyse Veffet d’un défaut de collimation de la lumiére et on présente les conditions exactes 
de convergence ou divergence de la lumiére. 


ZUSAMMENFASSUNG. — Der optische Wegunterschied zwischen zwei kohdrenten Strahlen bei der interferometrischen Registrierung 


einer Difusionszone wird bis zu einer Genauigkeit der dritten Ordnung fiir den Eintrittswinkel des Lichtes und fiir den gesamten 
Zuwachs des Brechungsindex in der Diffusionszone berechnet. Es ergibt sich, dass die Fehler dritter Ordnung einen Minimalwert 
fir diejenige Defokusierung der Zellenmitte aufweisen, die schon die Theorie zweiter Ordnung forderte, um die Wienersche 
schiefe Lage der Gradientenkurve zu beseitigen. Die Fehler dritter Ordnung sind stark vom Eintrittswinkel abhdngig und es 
gibt darunter einen bestimmten Winkel, der giinstiger als alle anderen ist. Die Fehler dritter Ordnung klingen mit der Zeit rasch 
ab und man kann eine gewisse kritische Diffusionszeit angeben, nach der die Fehler nicht mehr beobachtet werden kénnen. Diese 
kritische Zett ist fur den erwdhnten optimalen Eintrittswinkel ein Minimum. Die mathematischen Ausdriicke fiir diesen Eintritts- 
winkel und diese kritische Zeit enthalten nur bekannte Apparatkonstanten und die Zahl der Interferenzstreifen. Die Unschdrfe 
der Streifen, die von der Defokusierung der Zellenmitte und der vertikalen Ausdehnung der Lichtquelle herrithrt, wird ebenfalls 
mit der kritischen Diffustonszeit in Zusammenhang gebracht. Man kann sich gegen eine Auswertung von noch nicht fehlerfreien 
Interferogrammen schittzen. Das geschieht durch bestimmte Justierungen, durch die sowoht die Fehler als auch die Unschdrfe der 
mittleren Streifen gleichzeitig zum Verschwinden gebracht werden. Schliesslich wird auch der Einfluss von nicht parallelem 


Licht untersucht und es werden die mathematischen Bedingungen fiir die zuldssigen Betrdge der Konvergenz und D ivergenz des 
Lichtes aufgestellt. 


I. Introduction. — Interferometry has made _ it 
possible to measure diffusion coefficients with a for- 
merly unknown precision. It is carried out by taking 
photographs of interferograms from twin cells, one 
compartment of which is filled with a solvent through- 
out, while the other one contains the same solvent 
stratified on top of a solution of the substance to be 
measured in such a way that the starting boundary 
is very sharp. The interferograms taken during the 
course of the free diffusion between solution and 
solvent are interpreted by the simple equation: 


(1) AS=al[n(x)—ml, 


where AS denotes the difference in optical path 
length between the two coherent light beams, a the 


cell thickness, m the refractivity of the solvent, and 
n(x) the variable refracti'vity across the diffusion 
boundary in the cell (a is the vertical cell coordinate). 


As was pointed out in a previous article [1] on 
this subject, equation (1) is only a first approxima- 
tion to the real path difference. The equation would 
be exact if the light passed horizontally through the 
diffusion and reference cells, but it does not. The light 
may enter the cell under a finite entrance angle, and, 
if collimated horizontal light is used, it will in any 
case be deflected in the medium of variable refrac- 
tivity and describe a curved path, a trajectory, 
within the cell. A light pencil will thus sweep through 
a certain vertical region of the cell, within which the 
refractivity cannot possibly be regarded as constant, 
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especially with regard to the high accuracy obtain- 
able by interferometric means. As soon as the refrac- 
tivity along the trajectory varies more than corres- 
ponds to the accuracy of the interferometric method, 
the simple equation (1) can no longer be expected to 
be valid, and an exact path length difference can 
only be obtained by an integration over the length 
of the trajectory. There are also other reasons for 
deviations from the simple equation (1). For a more 
detailed discussion of this question, the reader is 
referred to the cited article. 

It was shown in the said article that the expression 
for the exact path length difference is too complicated 
to be of any practical use, but that it is possible, by 
a proper use of mathematic approximations, to 
derive equations which are practically useful and 
more exact than equation (1). Thus it was found, 
by taking care of the terms next biggest to the main 
term in (1), that the optical path length difference 
can be described by the equation : 


a8) (1 —2 7) n’ (2) 


Q) AS =a [n(2)—m] 4 5 + 
m 

a? (2 —3r) n' (x) 

6m : 


where r is a parameter open to free choice in the adjust- 
ment of the optical system. It was also shown that 
both the above aberration terms could be brought 
to disappear by certain tricks in the adjustment of 
the optical system. The last term, which is most 
serious and responsible for the skewness of the gradient 
curve predicted by Wrenur [2], can apparently be 
brought to disappear if r is chosen equal to 2/3. 
This corresponds physically to an optical adjust- 
ment making a plane 1/3 of the cell thickness from 
the exit wall optically conjugate to the photographic 
plate, 7. e. to defocusing the middle of the cell 1/6 
of the cell thickness. The first aberration term, on 
the other hand, can apparently be extinguished by 
making 5) = 0, but this is not necessary. Even in 
oblique, collimated, and in uncollimated light, this 
aberration can be cancelled by another optical 
trick, viz. by a certain coordinate transformation, 
the physical correspondence of which consists in the 
experimental determination of the photographic en- 
largement factor by placing the object in the plane 
of the middle of the cell, although this is not optically 
conjugate to the plate. More details of this, and of 
other interesting properties of the two aberration terms 
are to be found in reference 1. 

These results, which have been experimentally 
verified [3], are of course very favourable since they 
show that the simple equation (1) is worthy of a 
great deal more confidence than could be expected, if 
the tricks mentioned above are applied. On the 
other hand, they show that the second order theory 
is not sufficient to reveal the ultimate limit of inter- 
ferometry in the measurement of strong refractivity 
gradients. It was mentioned already in reference 1 
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that the second order theory could not give any infor- 
mation of the critical relation between the thickness 
of the cell and that of the diffusion boundary. There 
must be such a relation, beyond which a light pencil 
traverses too great a portion of the diffusion bound- 
ary, causing the optical path length to assume an 
average value over this portion rather than the desi- 
red path length along a horizontal line. To get this 
information, a third order approximation of the 
path length difference is necessary. It will be the 
object of this investigation to derive the third order 
aberrations, among which we will find the terms 
responsible for this averaging effect. 


II. Notation and basic equations. — The notation 
and the general scheme of the analysis will be 
explained with reference to figure 1, which shows 
the paths of two coherent light pencils through a 
diffusion cell and a reference cell. The z axis is the optic 
axis, and a monochromatic point source is situated 
at z = — p+ra, x = q. The distance p can then be 
regarded as a divergence parameter, p = « corres- 
ponding to collimated light, while the ratio g/p can 
be regarded as an off-axis parameter. Both parame- 
ters are open to free choice in an experimental arran- 
gement. The diffusion and reference cells are both 
situated between the planes z = 0 and z = a. The 
two coherent pencils are thus projected on to the 
same plane, although they may have any mutual 
orientation in reality. To begin with, the cell walls 
are imagined as infinitely thin. The diffusion cell has 
a variable refractivity, n(x), assumed to be inde- 
pendent of z, while the constant refractive index of 
the reference cell is called m. In order to facilitate 
the calculation, this medium is first assumed to 
extend even outside the cells, through the whole 
part of the optical system that we have to deal with. 
This assumption will be omitted at the end of the ana- 
lysis. The plane z = ra is assumed to be the Gauss- 
ian image plane of the plate where the interfero- 
gram is formed, thus the parameter r may be called 
the cell-focusing parameter. It is also open to free 
choice in an experimental arrangement. The optical 
imagery from this plane to the photographic plate is, 
moreover, assumed to be completely aberration-free. 
There is then no optical path length difference be- 
tween two light pencils passing from the same spot in 
the plane z = ra to the same point on the photograph- 
ic plate. This is the reason why our calculation of 
path length differences can be stopped at the plane 
Sh VIT Os 

One coherent pencil enters the diffusion cell at 
xz == 2%» under the entrance angle 8, describes a 
curved path within the cell, leaves it at « = x, under 
the exit angle 8,, and strikes the plane z = ra at the 
level « = w,,. The internal entrance and exit angles 
are called «, and «,, respectively. The other coherent 
light pencil takes, due to our simplifying assumption, 
the direct rectilinear course from the light source to 
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Fic. 1. — Light paths from light source to the left through reference and diffusion cells to the Gaussian 
image plane, z = ra, of the photographic plate. The cells are situated between the planes z = 0 and z = a. 
The light path through the reference cell is rectilinear, whereas that through the diffusion boundary is curved. 
The two pencils are coherent since they originate at the same spot of the light source. Moreover, they are ca- 
pable of interference on the photographic plate, since they both arrive at the same spot of its optically con- 
jugate image plane. The plane z = a/2 is used in the calibration of the cell(in the determination of the pho- 
tographic enlargement). Therefore, the coordinate x,,, and not 2,q, is the coordinate read on the photogra- 


phic plate. 


the same point w,, in the plane z = ra. This pencil 
strikes the plane through the middle of the cell 
at © = «&,. 

It is now easy to formulate an exact expression 
for the optical path length difference between the 
two coherent rays. It is 


— "a n (x) dz 
ae mee a (2) 
cos 38, o COS « 
__ ma (7—1) mp 
' cos 8, cos 8,° 


As already pointed out, this expression is of no prac- 
tical use since it contains an integral of a function 
that is not directly integrable, and since it also 
contains the angle 8,, which is not known to an 
experimental worker. The work now to be performed 
consists in the computation of a sufficiently accurate 
approximate value of the integral, and in the elimi- 
nation of 8, and other unknown variables that may 
appear in the integral. Moreover, the whole expres- 
sion has to be converted into a convenient form, 
suitable for drawing practical conclusions regarding 
the best choice of the free parameters. 

Each term in equation (3) will be submitted to 
mathematical approximations, and since this article 
is called a third-order theory, we will then take into 
account all terms being one or two orders of magni- 
tude smaller than the main term given by (1). The 
concept of order of magnitude then apparently 
requires some further explanation, 


If equations /1) and (2) are divided through by a, 
we get equations for the relative difference in optical 
path length, and both sides of the equations become 
dimension-less. The first approximation of this rela- 
tive difference is identical with the refractivity differ- 
ence, n(x)—™m, at the levelx. As long as we are concern- 
ed with dilute solutions, and as long as we are using 
the solvent in the reference cell, this is a quantity 
that is small compared to unity. In general, it is not 
greater than 0.002, but we may allow it to rise to the 
order of magnitude of 0.01. By definition, we will regard 
a dimension-less quantity of this order of magnitude 
as a quantity of the first order or first degree. The 
total refractive index increase across the diffusion 
boundary will be denoted An and will be assumed to 
be of the same order of magnitude, 7. e. also as a quan- 
tity of the first degree. 

All angles, 8, 8,, 8,, %, and «,, will be regarded 


as first-degree quantities on this scale, i. e. they 
will be assumed not to be appreciably larger than 
0.01 radians. 

The quantity an'(z) is identical with the angular 
deflection of a light pencil in the refractive index 
gradient, and will thus also be regarded as a first 
degree quantity. Now, the first derivative can be 
written, for an ideal diffusion boundary : 


An 


4 '@) =" fy), y= ap, 


where f(y) is the Gaussian error function, and bp is 
half the distance between the two inflexion points 
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of the normal error curve. The numerical values of 
the error function lie between 0 and about 0.4 and 
must consequently be regarded as quantities of the 
0: th degree. Since An has already been defined as a 
first-degree quantity, it can be concluded that a and 
b must be of the same order of magnitude. This is a 
reasonable conclusion, standing in conformity with 
general experimental practice. 


The angular off-axis position of the light source, 
q/p, will, as all other angles, be regarded as a first- 
degree quantity. The same assumption is made con- 
cerning the angular extension of the diffusion gra- 
dient as seen from the light source. Thus 2/p is also 
a first-degree quantity, as well as b/p. Further, since 
a and 6 are of the same order of magnitude, the quan- 
tity a/p, if it appears, will be treated as a first-degree 
quantity, although it cannot be physically interpret- 
ed as an angle. 


Finally, the second derivative, n” (x), will appear in 
our approximations. Its order of magnitude must 
be the same as that of the first derivative, since it 
contains An as a factor, and since the first derivative 
of the error function is of the same order of magni- 
tude as this function itself. It is to be noted, however, 
that the derivatives of n with respect to x are not 
dimension-less. They bave to be multiplied by a 
distance in the same power as the order of the deri- 
vative under consideration before they can on the 
whole be compared to other quantities. The distance 
that appears in combination with the derivatives is 
the cell thickness, which, as already agreed upon, 
is regarded as being of the same order of magnitude 
as the thickness of the diffusion boundary. 

To sum up, we can thus state that, on simplifying the 
different terms in equation (3) by mathematical approx- 
imations, we have to take care of all terms up to 
and including terms of the third degree in any com- 
bination of the above-mentioned variables of the 
first degree. In order to examplify this principle, we 
can start by inspecting the two aberration terms in 
equation (2). Since r is a dimension-less parameter 
of unity order of magnitude (0: th degree), both these 
terms are found to be of the second degree (after 
division by a), since the first one contains 8) and 
an'(x) and the second one a?n'*(x). A 3, term would 
have to be included in (2) if it existed, but the 
mathematical analysis in the previous article showed 
that such an aberration does not exist. The fact that 
the two aberrations in (2) are of the second degree 
referring to the quantitative scale agreed upon here 
is of course the reason why the corresponding errors 
were called second-order aberrations. 

In the third-order theory to be presented here, 
we thus have to take into account such terms as 
a(n—m) &, a 3 n(x), a5 n'* (x) n"(x), at n'* (2)/p, 
to mention a few examples of terms that were neglected 
in the second-order theory. 


The different variables are interconnected by the 
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following equations, the validity of which is easily 
realized. 


(5) fon se ee 
p—ra 

(6) Mm 8in 89 = n (Xo) SiN a, 

(7) n (%q) cos «, = n (X_) COS &%, 

(8) n (z,) sin «, = m sin 8,, 

(9) Lq = ©, + a (r — 1) tan 3,. 


By squaring and adding equations (6) — (8), we 
get 
(10) m? (sin® 5, — sin?’6,) =n? (@,) 1" (a). 


Along the trajectory, we have the following equa- 
tions : 


(11) dS Bg (x) 
dz COS a 
(412) n (x) cos « = a constant, 
(13) a == Ube 6 
da n' (x) 
4 — 
at dz n (x) 


Equation (14) is obtained by differentiation of (12) 
with respect to z. 


III. Second-order approximations of auxiliary 
variables. — In order to carry out the computation 
of the third-order approximation of the path length 
difference, it is necessary first to derive the second- 
order approximations of some auxiliary variables. 
Thus, the only vertical coordinate that is known to 
the experimental worker is that on the photographic 
plate, and consequently all other coordinates such as 
x, and x,, have to be eliminated from our equations 
in order to make them practically useful. The same 
thing can be said about the exit angles, 8, and «,, 
as well as the internal entrance angle %. 


The exit coordinate, x,. — In order to derive an 
approximate expression for this coordinate, we will 
so far regard x) as known, since the only thing that 
can be done at the start of the analysis is to estab- 
lish a relation between the entrance and exit coor- 
dinates. Following the curved light path through 
the cell, we can regard x as a continuous function of 
z in the interval between z = 0 and z = a and apply 
a MacLauvrin series: 


1 Oe a /a*a\ 
15) 2,=2,+°2/(% a a Ete 
ce va Of (ce Maes re a 
iar i Ri) 2 a* /d4z =a a® /d®x 
Gece, 24 i mii? 120 cae 


where the last term is the residual term with an un- 
known numerical figure 6 lying between 0 and 1. It 
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is to be noted that the relation (15) holds whether 
the power series in a converges or not. Thus we 
need not worry about any convergence criteria for 
the power series, if only we can show that the resi- 
dual term is of insignificant magnitude on our scale. 
The first derivative in (15) is taken directly from (13), 
and the successive differentiation is easily carried 
out with the aid of (14). Thus we get the following 
second derivative 
d?x hat 
ei aa a) 


(16) 


where n’ and n are abbreviated symbols for n‘(2) 
and n(x), respectively. This abbreviated notation 
for n and its derivatives with respect to 2 will be 
used consistently in the following. In the above 
equation tan « as well as n’ (together with a, which 
appears as a factor in (15)) are of the first degree on 
our scale of orders of magnitude. Hence the second 
term, being of the third degree, can be omitted since 
we are looking for only the second-degree approxima- 
tion of x,. Thus, we have 


Ge 
(17) a ae 
Equation (17) can also be used in the computation 
of the higher derivatives, because the third-degree 
term in (16) does not give rise to any terms of lower 
degree on further differentiation. 

The next differentiation of (17) gives, with the aid 
of (13), 


3 Wt L2 
(18) a =(4 =! ) tan a, 


n n? 


without any approximations. Here again, we can 
discard the second term, since it is of the third degree, 
and we get 


19 ee o n! 
eo dz? n 


because terms of the third and higher degrees in tan 
« can also be disregarded. Differentiation of (19) 
once more gives only one term of the second degree 
déx et 

20 = 
eu dz* n* 


The fifth derivative contains, as is easily realized, 
only terms of the third and higher degrees. The 
residual term in (15) is thus too small to be considered 
and can be omitted. Consequently we get the follow- 
ing second order approximation of the exit coor- 
dinate : 


a? ni aan. 

(21) eer. np ge ete vee 
‘2 g Z DNs 6 no 

4 t MW 

a* n,n, 

2 

24 ny 
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The external entrance angle is known to an expe- 
rimental worker, whereas the internal one is not. 
Thus we will eliminate « in favour of 8) with the aid 
of the following relation, easily realized from equa- 
tion (6) : 


MS8o 
(22) Coney 


+ correction terms of the 


No 


third and higher degrees. 


The correction terms can of course be omitted in a 
second order approximation, and we get 


Sel 3 " 
23) mMa8, ee Ma* SoM 
“L, = «£, 

( a or al No Qn Gn 
0 

4 Ly, 

a*nyn, 

2 

24n, 


It is, however, inconvenient to have the variable No 
in the denominators. Therefore, it is advisable to 
make use of the following relation, 


(24) = —(- =e m 

Ng m m No m MNo 
ee (1 fey me), 

m No 
Hence it is possible to replace 1/n, by 1/m if a correc- 
tion term of the next degree is added. (According to 
our assumptions regarding orders of magnitude, no-m 
is a quantity of the first degree). The same relation 
can, of course, be applied over and over again to the 
correction terms, which still contain n> in the deno- 
minator. If, now, relation (24) is applied to the two 
first-degree terms in (23), two additional, 7. e. in all 
four, second-degree terms are obtained. All of them 
also contain m» in the denominators. The relation (24) 
can be applied to these terms also, which results in 
the formation of four terms of the third degree. We 
are, however, dealing with a second order approxi- 
mation only, so the last-mentioned terms are omitted. 
Consequently, we are allowed simply to replace 1/n, 


by 1/m in all terms of the second degree, and thus 
get the equation : 


an, a8, (n»—m) 
(25). &, = @, ++ a8, + Sc el = 
2m m 
a a’n! (n,—m) n aS Ng atngn, 
2m? 6m 24m? 
The internal exit angle, x,. — This angle will be 


derived in the same way, thus we put 


ee ae Heal 
“a “or fl (Ss), ¥ 2 (S)-. at 5 (as). 


and start the successive differentiation from equa- 
tion (14), which gives 
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d2a an” da nin" 
i ; =e a 

dz? n dz? n? 


Higher derivatives only contain terms of the third 
and higher degrees. Thus we get 


an, aza,n" a®nin’ 
(2.8) et mae P| z J nee ae “= 
No ano 6n, 
The external exit angle, 8,. — It should be pointed 


out here that this angle is external with reference to 
figure 1, where r> 1, but internal for r <1. Its value is, 
however, independent of r, and is calculated from 
equation (8). Since we are dealing with a second 
order approximation only, we can put the sines of 
this equation equal to their angles, and, further, by 
the use of a first order TayLor expansion of n (2,), 
we get 


(27) m3, =a, [2 (%) + (%g—29)n' (25) |- 


The factor «, is of the first degree. According to 
equation (25), that is also the case with («,—w,) after 


division by a, and an'(a 9) is, according to our defini- 
tions, a first-degree quantity. Consequently, the 
second term in the above expression becomes of the 
third degree and has to be disregarded in our second 
order approximation. We thus arrive at the simple 
relation 


(28) ee 


me 


With the aid of (26), (22), and (24), we get the follow- 
ing second order approximation for the external 
exit angle: 


/ " 
an ar?on 
0 o°"o 
= ue 4 | 
m 2m 


Sanat 
a®ngno 


6m? 


(29) 8 


The coordinate in the Gaussian image plane of the 
plate, x,,. — This coordinate is given by equation (9), 


in which tan 8, is approximated by 8 


ra‘ 
» and in which 


we then introduce the expressions for x, and 8, Just 
derived. The result is 


WNg (2r—1) a8o (NMp—mM) 


Oa do am 7 = 
a?n ((No—m) a®3,n, (3r—2) atngn,(4r—8) 
—"s | 
2m? ; 6m 24m? 


The coordinate of the reference pencil in the middle 
of the cell. We have a considerable interest in this 
coordinate because of the conclusions drawn in the 
second order theory in reference 1. We will start this 
calculation from the exact equation: 


} Lig 4 
9 3 2 
(31) a (r—1/2) P 
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By adding and subtracting xz) in both numerators, 
and with the aid of (5), this equation can be written 
in the form 


(32) 2,,—%) = XLq—Xy — ee — 
a(2r—1) (2, —1a,) 
ae. xp 


‘The second order approximation of the first two terms 
above is given by (30). The third term is of the first 
degree and is added to the corresponding term in (30). 
The last term in (32) is of the second degree, since 
a/p is a first-degree quantity. In this term, it is conse- 
quently sufficient to introduce the first-degree terms 
of (30) exclusively. These substitutions give the 
result : 


as, an, (2r—1) a, (ny —m) 


(33) i |. - — 
u Y « T 
ae 2 2m m 
ee ; " : 
Gligiitp mt) a®3n, (3r—2) 
¢ | > 
2m? 6m 
if 
asnin’, (4r—3) a®n} (2r—1)? 
24m? Amp 


IV. The third order approximation of the path length 
difference. — We are now ready for attacking the 
different terms in the exact equation (3), and we will 
start this work by studying the path length within 
the cell. 


The optical path length of the trajectory. — This 
length is easily obtained by a MacLaurin approxi- 
mation of the same kind as used earlier. We can write : 


2 a n(x) dz a /as 
(34) te bens | NK ) ec (=) xi 
JO CON { \dz/z=0 
ae ee ei ae a i at (Sa) L 
F Ndzi7 a2 6.0 6) \dz2) g0p en DEG Nae eee 


ene ee 
ee ary ie ot 


and omit the residual term A, if it can be shown to 
be of insignificant order of magnitude on our scale. 
The successive derivatives are obtained by differen- 
tiation of equation (11) with respect to z, making 
use of the relations (13) and (14). In each derivative, 
we are allowed to omit terms of the fourth and higher 
degrees. We then get 


ds d2S d3S on* 
yn (! 2/9, — f = 2 «2 " 
dz es: dz? a dz® n pale 
(35) 


d4s San'n" d5S 8n'2n" 
dz ith » toed n2 


The residual term contains the sixth derivative of S 
with respect to z, which is of the fourth degree in An 
and «, and need not be considered. Thus we get the 


Q 
Vv 
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following third order approximation of the optical 
path length of the trajectory : 


2a n(x) dz mas, , 
(36) 8S, A De ang 4 Tee. 
ma? 3M an m®a®8,n, mats,nin, a ony th 
= oF came. - 2M ae 
No 3No oN, 3N, 15n, 


Before we proceed to the next term in equation (3), 
we will, however, eliminate the variable mn) from 
all denominators by application of equation (24) to 
three terms of the second degree, and by putting 
1/no = 1/m in all terms of the third degree. We then 
get an equation with six terms of the third degree, 
which, however, is more convenient to use in the 
following treatment : 


ra n(x) dz mas, 
om ae : ae a 
(ii = | ae dk, ao ian 
Jo cose 2 
12 Q2 
arn, a8)(ng—m) 
+ a®s,n, 4 
3m 2 
2 ! ome. 3,,/2 mF 392, 
a*3,n(n-—m). an, (nom) — BSN, 
= = 2 = r 
m om? 3 
a le Geshe med 
: L* 8 Nolo A°Ny Ng 
3m 15m? 


The optical path length from the exit wall to the 
Gaussian image plane of the plate. This length is 
represented by the third term in equation (3), and is 
easily calculated, since we have already the second 
order approximation of the angle 8,. We have 


38 ma (r—1) = 
Pe cos 3, 


ma (r—1) (1-+82/2) 
in which equation we have to introduce (29). If the 
latter equation is written in the abbreviated form 


(39) p= 8a + 8459 

where the first term represents the first-degree part, 
and the second one the second-degree part of 8,, then 
the third-degree approximation of its square can be 
written : 

(40) soo 288 

Thus the computation of (38) is fairly easy, and the 
result is 


ma (r—1) mas, (r—1) 
4 Bee Ss , | 
(41) Pa mra — ma 4 5 4 
’ atn.” (rj) 882" (r—1) 
+a?8,n, (r—1) + a5 : 3 S. 
2a*dngng (r—1) ann’ (r—1) 
3m 6m? 


H. SVENSSON 


[Opr. Acta 


The optical path length from the light source to the 
entrance wall, and the total path length of the reference 
pencil. — These quantities, which are represented 
by the first and last terms of equation (3), respectively, 
are preferably treated together. We have 


m (p—ra) mp 
42 S,— 84 = - 
oe : 5 COS So cos}, 
mp (cos8, — cos8p) mra 
COSS cosd, cosd) | 


According to figure 1, we have, however, the rela- 
tion : 


(43) Lra—] = Xpg—Xq pig g = Urq Xo iis 
+ (p—ra) tans, = p tans, 
from which we get, without any approximation, 


to, Fe tana, 


4A — 
Ces P tand, — tandy 


Introduction of this value of p into (42) gives 


(45) 


SS ne re 
cosd, — coss, mra 
—ra tan8,) Be: : 
sin (8,—55) COS 8p 


We know from equation (30) that x,,—a, is of the first 


degree and conclude that the trigonometric function 
of the two entrance angles has to be computed to 
the second degree. Hence the first-degree approxi- 
mation of the sine function and the second-degree 


approximation of the cosine function are adequate, 
which gives 


(46) S,—S, = =m — 2, —1as,) (oa 


— mra (1+ 85/2). 


Since the trigonometric function also appears to be 
of the first degree, it suffices to use the first-degree 
approximation of the tangent function in the first 
factor. It now remains to express the entrance angle 
of the reference pencil in terms of that of the pencil 
passing through the gradient. The accuracy required 


is that of the second degree and the desired relation 
becomes: 


Nb] Comal Lg Xo 
OF — =— — 


P. p P ie 


ra 


+ (1—ra/p) tan 8 = 8, -4 


0 


P 


Introduction of the first-degree terms of (30) further 
gives 


a*n, (2r—t) 


La 
(47) ae 


8 = % + 


After introduction of this value into (46), the latter 
equation can be written : 
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ani, (2r—1) 
/, ag Alas! by PO ae 0 7 ” 
(48) S,—S,= m am =. (ute 
F a*n’, (2r—1) mrad? 
8 + — mr a ., 
a 4 mp 


where the subscript 2 is used to denote the second- 
degree terms of the function. On multiplying the two 
parentheses together, it is observed that the pro- 
duct of the two second-degree terms can be omitted, 
and the result is 


52 2 ! 
(49) S,—5.= mrasy a 8514 (2r—t ) 
ey ig —— NG mra— —, = “u 
; a?8n5 (n>—m) a®8ini, (3r—2) 
+ a8, (n,—m) -— ———_——_ - - 
2m 6 
Lt [pat 4y!2 (D» 2 
a*8 nn, (4r—3)-  atn, (2r—t) 
24m 8mp 


The total path length difference. — This is now 
easily obtained by adding equations (37), (41), and 
(49). The result is 


a*8.n an.” (3r—1 
(50) AS = a (n——m) + es Pesce es 


a8.(ng—m) a8, (Ny —m) 


0 
eo 2m : 
ain!” (ng—m) aon), 
ePesoilae. aio 
a8,n,n, (12r—d) ann, (5r—3) 
"Ban Bm 
atn,” (2r—1)? 
Smp 


This equation would be immediately useful for our 
purpose if z were directly related to the vertical 
coordinate on the plate carrying the interferogram, 
but it is not. The reader is reminded about the fact 
that the symbols n,, nj, n, are abbreviated forms 


of n(%), n'(%), and n(x), respectively, and that 
the entrance coordinate x, is not known to an expe- 
rimental worker. No conclusions of interest can 
therefore be drawn from equation (50). 


VY. The coordinate transformation. — Physical inter- 
pretation of coordinate transformation. — After it had 
been shown in the second order theory that the 
skewness aberration could be brought to disappear 
by defocusing the middle of the cell 1/6 of its thick- 
ness in the direction towards the light source (by 
making r = 2/3), it was also found that the remain- 
ing second-order aberration, in the case of a light 
source without vertical extension, could be cancelled 
by experimental determination of the photographic 
enlargement factor using an object placed in the midd- 
le of cell in spite of the defocusing mentioned. When 
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this is done, the cell should be filled with a medium of 
the refractive index m, and the light source system 
has to be the same as that used in the actual diffusion 
experiments (the same divergence parameter p, and 
the same off-axis position parameter, g/p). In prac- 
tice, the diffusion cell is not at all used in this opera- 
tion, but a transparent glass scale is placed on the 
cell stand at the same optical distance from the lens 
system as the middle of the cell in the actual experi- 
ments (correction for the cell wall and for half the 
thickness of the reference cell if the scale is not sur- 
rounded by a medium of refractivity m). The same 
procedure applies if the optical system is not free 
from distortion, that is, if the photographic magnifica- 
tion is not constant. The only difference is that one 
should talk about a calibration of the plate coordi- 
nate instead of a determination of an enlargement 
factor. 

It is easily realized with reference to figure 1 that 
the procedure recommended makes the coordinate 
L,,, the coordinate of the reference pencil in the middle 


of the cell, directly related to the plate coordinate, 
since no refractive index gradient is used when the 
calibration is carried out. We thus conclude that it 
is impossible to interpret equation (50) without eli- 
mination of the coordinate «, in favour of z,,. 

If we denote, by definition, the right-hand side of 
equation (33) by Az, we thus have to introduce 

(51) L, = t, — Ax 
in all terms of equation (50) where a» appears. This 
is the reason why it was advisable to get rid of n(2%p) 
in all denominators. However, we now encounter 
another difficulty in the fact that Aa itself contains 
the variable z) which it shall serve to remove. Thus 
we have to start with removal of 2, from the right- 
hand side of equation (33). 


Removal of the entrance coordinate from the equa- 
tion for the plate coordinate. — Equation (33) contains 
two terms of the first degree and five terms of the 
second degree on the right-hand side. If we take an 
arbitrary term of the last kind, e. g. the last term in 
equation (33), and introduce the new independent 
variable with the aid of (51), we get 


a®n, (2r—1)? a? (2r—1)? n'(”,—Ax) = 
aie Amp m 
eat — he 


ju! (%,,)-—Ax ney) | 


Amp 


The second term in the parenthesis is one order of 
magnitude smaller than the first. Since already the 
first term is of the second degree, the second one must 
be of the third degree and need not be included into 
a second order approximation. We conclude that 2% 
can simply be replaced by z,, in allsecond-degree terms 
of (33), without impairing the exactitude of the equa- 
tion. By the same reasoning, it can be concluded 
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that it is sufficient to apply the first order approxi- 
mation of (33) : 
a8 a? n, (2r—t) 


i Ah = 


} x 
0 m y) 2m 


in the removal of x) from the two first-degree terms 
of the same equation. Now, the first of them does 
not contain Zp, so the substitution becomes in fact 
restricted to one single term : 
a’n, (2r—1) a? (2r—1 
(53) : = ( ) n'(x,,) — 


2m 2m 


“m 
a*n,, (2r—t) 
2m 


atn’ nn” (2r—i)? 


mom 


8 "( 
- n(x 
2) m 


" 
a*Son,,, 


nl (2) 


(2r—1) (2r—1) 


aan’ 
2m 4m 


“Mm 


4m? 


Heren}, andn’, are abbreviated forms of n‘(a,,) and 
n"(x,,). If (53) is introduced into (33), and if 2, is 
changed to z,, in all second-degree terms, we get the 
following equation to be used in our coordinate 
transformation : 


(54) £, = &,, — (Ax), — (Az), 


m 
where (Az), and (Az), represent the first-degree and 
second-degree terms, respectively, and are given by 
the equations: 


a8p an}, (2r—1) 
A — 
as, (n,,—m) an’, (n,,—m) 
6 A an 0 m m m 
(56) (Ax). om Som 
a®%8,n), ain! n), (24r2—28r+9) — a8n/ (2r—1)? 
12m 24m? ~— 4mp 


Removal of the entrance coordinate from the expres- 
sion for the total path difference. — We are now ready 
to attack equation (50) with the coordinate trans- 
formation (54). For reasons quite analogous to those 
presented in the foregoing section, it is superfluous to 
transform terms of the third degree; in them, 2 is 
simply replaced by w,,. In terms of the second degree, 
it is sufficient to use the simple relation : 

(57) n'(X%,) = n'(x,,) — (Ax), n"(2,,) 
whereas in the only existing term of the first degree, 
we have to use the more complicated equation : 


(58) n (%) = n (2) — |) shi | 1" (Bq) 
(Ax); 
+ 9 n"(x,,)- 
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The substitution involves no special difficulties, and 
one gets the following result: 


3, | ad,n' ml 
(59) AS = a(n—m) | 1 + 5 oe ar 6m? 
a®n'® (2—3r) 7 eon" as nH 
6m Ps 24m 
_ an'n"(15r3—20r-+7) atn'?(2r—1)? 
. 30m? 8mp 


The effect of plano-parallel walls. — At the start of 
this analysis, two facilitating assumptions were made, 
first the assumption of infinitely thin cell walls, and 
second the assumption of an external refractivity 
equal to that of the reference cell. The first assump- 
tion is equivalent to the removal of two plano-parallel 
plates of a thickness equal to that of the walls, where- 
as the second one corresponds to the insertion of a 
plano-parallel plate of thickness p—ra between the 
light source and the diffusion cell. These assump- 
tions, which are thus of the same nature, facilitated 
the calculation greatly, since thereby the path of the 
reference pencil became rectilinear. 

We are now going to free ourselves from these 
unrealistic assumptions simply by replacing 8p, 
whereever it appears in (59), by 8)/m, by replacing 
p by mp, and by redefining the distance p as the optical 
distance between the light source and the Gaussian 
image plane of the plate. This is the procedure gene- 
rally adopted in elementary optics, and leads to the 


equation : 
3, as,n' ene 
== es 
(60) AS=a(n—m) | 14 Poa a7 anor ae 
| an’? (2—3r) — a8n" atsn'n" 
| 6m 2hm? ' 24m? 
ann" (15r2—20r+7)  atn’? (2r—4)? 


30m? 8m*p 

The correctness of the above procedure is of course 
not self-evident in a third-order approximation theory. 
It will be shown in Appendix 2, however, that the 
errors involved are negligible. 

Equation (60) is thus our final equation for the 
optical path length difference. As is readily seen, it 
contains only one aberration term of the second 
degree. In the second order theory, two such terms 
were found as long as the independent variable was 
the coordinate x, in the Gaussian image plane of the 
plate. The term containing a8 )n' was then shown to 
disappear by a second coordinate transformation. 
In this treatment, we have transformed directly to 
the coordinate «,, of the reference pencil in the middle 
of the cell, so the absence of an a8,n’ term in (60) 
actually proves once more the correctness of the con- 
clusions drawn in the previous paper. Since in the 
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following treatment we have no other vertical coordi- 
nate to deal with than z,, the subscript has been 
dropped in equation (60). 


The optimum position of the cell. — The obvious 
conclusion that the only remaining second-degree 
aberration disappears for r = 2/3 was drawn already 
in the previous article. There are only two more 
terms that depend on the longitudinal position of the 
cell, and among them the last one is present only 
in uncollimated light, which is comparatively seldom 
used. It is important to observe, however, that 
the a> term acquires a minimum value exactly for 
r = 2/3. Consequently the third order theory gives 
additional evidence in favour of the cell position 
recommended in the second order theory. We will 
consequently use this particular value of the parameter 
r in the following analysis, and equation (60) then 
takes the form: 


; . | de 0 OS an! 
(61) AS = a [n (x)—m] | 1 4 ami = aa + 
oe eo as nin ann! asn'? 
Gat > tee + tome boone + Famip 


VI. Physical interpretation of the third order approxi- 
mation. — If equation (61) is written in the form 


(62) AS = a [n(x)—m ] [1 a R(a, t, So) | - A(z, t, 30), 


it is evident that we have to do with two different 
aberration functions, one absolute, A (a, t, 39), and 
one relative, R (x, t, 8). It is advantageous to treat 
these two functions separately. 


1. The absolute aberration function in collimated 
light. — In collimated light, we have an infinitely 
large p value, and 8) is independent of x. The abso- 
lute aberration function reduces to three terms : 


aen'?n" 
90m?r 


Here we have divided through by \ in order to get 
the aberration in terms of light waves. 


ae syn" atds,n'n" 
A(x, t, 89) = + : 
m4 24mr 24m*r 


(63) 


Light source on the optical axis. In this case 8 = 0, 
and the aberration function reduces to one single 
term. It is of great interest to investigate the magni- 
tude of this term, since horizontal, collimated light 
is frequently used in practice. 

The refractive index function prevailing in the cell 
has already been assumed to be that given by an 
ideally diffusing substance: 

a er hy) 

2\/ xDt b 
where / (y) is the error function and 0} half the dis- 
tance between the two inflexion points of the curve 
n(x) ; 


(64) n’ (x) 
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The quantity An will be regarded as positive. Hence 
the error function f(y) has the same sign as n(x), 
which is negative in most practical applications. It 
follows further, since z = by and dx = 5b dy, that the 
second derivative is given by the equation 


(65) 


An An 


jal a 


The last term in (63) now takes the form 

CRM GEE et et 

90m? br 
This quantity is at a maximum for y = + 1/3, in 
which point f?f’ has the value 0.022235. The precision 
that can be reached in interferometry without using 
photometry is of the order of 2/50. Consequently we 
must require that the maximum of the above expres- 
sion be less than 0.02, which leads to the condition 


b* > 0.01235 (4m)? 


(66) sna) = 


(67) A (a,t,0) =a 


mn 


(68) 


If we now introduce the total number of fringes in the 
interferogram 


(69) Ni aan, 
the condition can be put into the form 
ra V3? \/ rat 
(70) 5? > 0.1111 ey b > 0.3334 ——=—. 


m 


Since 5? is proportional to the time, this condition can 
be regarded as a critical time condition, the aberra- 
tion (67) being smaller than 1/50 of a light wave after 
the elapse of a time defined by (70). 


Light source in an off-axis position. In this case, the 
entrance angle has a finite, but constant value. In 
order to facilitate the mathematical analysis of the 
aberration (63) as much as possible, the variable y 1s 
first introduced instead of x with the aid of the rela- 
tions (64) and (66). Two new, dimension-less para- 
meters B and C defined by the equations : 


2 45: Vane 


Vn | 
ymanrlt 
(72) 8 = C —— 
Va 
are then used for elimination of 6 and 8). The expo- 
nents to the right in these equations have been chosen 
in accordance with equation (70) what concerns 3, 
and for 8) in such a manner that the parameters m, 4, 
a, and N disappear altogether. Equation (63) 
consequently takes the form 


19 / 
pees =; 
= apt 


(71) b 


Cit ome 


(73) A (Y, B, C) “OAR a \D: 
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We will now investigate the conditions under which 
this aberration is smaller than 1/50 of a wave over the 
entire interferogram. 


The aberration as a function of the entrance angle. 
The partial derivative of A with respect to C is 


oA 2Cf' hs 
ac 24B®  * 24B3 
This derative vanishes for 
(75) C= — 7/28 
or, if we go back to the original variables, for 


(76) So — an'/2. 

Since an‘is the angular deflection in the cell, the con- 
clusion can be drawn that this aberration is at a mini- 
mum when entrance and exit angles are equal but 
for the sign. This fact was revealed by Loncswortu [4] 
several years ago by experimental tests with the 
Schlieren-scanning arrangement. It is also easily 
conceivable and plausible, since the w interval in the 
cell swept through by the ray must be at a minimum 
for a symmetrically passing ray. It is also interesting 
in this connection to compare the refractive index 
gradient in a plano-parallel cell with a prism of const- 
ant refractive index. The former gives a maximum 
deflection, whereas the Jatter gives a minimum deflec- 
tion to the symmetrically passing ray. 

It is also of interest to see how astonishingly small 
this minimum aberration is. If the value of C accord- 
ing to (75) is introduced into (73), the function A 
is found to assume the value 
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This aberration is 16 times smaller than the corres- 
ponding aberration for horizontally entering light, 
given by (67). However, we have not very much use 
of this minimum aberration, because experimentally 
it is not very easy to make the entrance angle a func- 
tion of the level in the cell such as to satisfy condi- 
tion (76) in every point. Moreover, we have the time 
dependence of n’ to take into a ccount. 


(77) =. 


** min. 


The aberration as «a function of tume. Next we will 
consider the partial derivative with respect to B: 
(78) aA Cae Chf Zi te 
fe . See - = — 
oB 12B° 8B4 45B° 


if we put this derivative = 0, an equation of the 
second degree in B is obtained, the roots of which 
are 


(79) 


By = 05702 9/0+B = = 20sec 


Negative B values lack physical interpretation. 
Consequently the roots (79) have no physical mean- 
ing if C and f have the same sign. The function A 
is then monotonically decreasing when B increases 
from 0 to infinity. On the other hand, if C and f have 
opposite signs, the function first decreases to a mini- 
mum, then increases to a maximum, and finally falls 
off to 0. The values of this minimum and this maxi- 
mum are 


(80) ATE = — 0.008493 C* y/f, 
A. = — 0.011229 C+y/f. 
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Fic, 2. — The absolute aberration A in collimated light as a function of time (B 5 i 
J a : at some differe 
parts of the interferogram, and for BC = — 0.1195. The values x = b and y ae ate to “ihe 


inflexion point of the gradient curve. Each curve has a minimum and a maximum. In the central parts of 


the interferogram (y < 1), these extrema occur comparatively late and are not ve 
parts of the interferogram (y > 1), they take slags at an : MEN Mircea 


c e the outer 
early stage of the diffusion process, and the minimum, 


especially, is very sharp. After the critical time, represented by B = 0.21, the aberration is below the 


resolving power for all parts of the interferogram. 
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Consequently there is a constant ratio between 
these two extrema (1.322). The ratio is rather close 
to unity. This implies that the aberration, for fixed 
values of C and y, during its rapid decrease in the 
course of time, stays for a rather long period at the 
order of magnitude given by (80). It is interesting to 
note that the increase of the aberration from its mini- 
mum to its maximum occurs shortly after the mini- 
mum aberration corresponding to a variable entrance 
angle has been reached [compare equations (75) and 
(79)]. The decline of the absolute aberration in the 
course of time can be studied in figure 2 for some 
different parts of the interferogram. 


The aberration as a function of the level in the cell. 
The time parameter B and the entrance angle parame- 
ter C will now be considered as constant, and y as 
variable. This function of y cannot be treated ana- 
lytically, but instead we can study the family of 
curves presented in figure 3, showing the aberration 
function through half the diffusion boundary for some 
different values of the product BC. It is seen here 
that the aberration is always zero at the centre of 
the gradient, then rises to a maximum, falls to a 
minimum, rises to a second maximum, and finally 
falls off asymptotically to zero for large y values. For 
small negative values of BC, the first maximum is 
large and the second one small, and vice versa for 
large negative values of BC. Consequently, when it 
is required to find the most favourable conditions 
with regard to the whole interferogram, the two 
maxima should have the same value. This occurs 
for 


(81) BC = — 0.1195 


Fic. 3. — The absolute aberration A in collimated 
light as a function of the cell coordinate. The centre 
of the diffusion boundary is at y = 0 and the inflexion 
point of the gradient curve at y = 1. The function 
has, in each half portion of the interferogram, two 
maxima and one minimum. For a special relation 
between entrance angle and time, given by the 


equation BC = — 0.1195, the two maxima have 
the same value. Using the critical time shown in 
figure 2, one thus arrives at an optimum entrance 


angle. 
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The absolute aberration A in collimated 
light as a function of the cell coordinate. The centre 
of the diffusion boundary is at y = 0, and the 
inflexion point of the gradient curve at y = 1. The 
figure shows the aberration at the optimum entrance 


Fic. 4. 


angle, BC = — 0.1195, in comparison with the 
aberrations for horizontally entering light, BC = 0, 
and for the optimum entrance angle of the wrong 
sign, BG = + 0.1195. 


and the two equally large maxima then have the 
value 


(82) Amax, = 0.01421/360 B4, 

The curves in figure 3 all belong to negative values 
of BC. In order to give an idea of the magnitude of 
the aberration for non-negative values of BC, figure 4 
has been constructed. Here the aberration curve 
for BC = 0 (horizontally entering light) can be com- 
pared to that for BC = — 0.1195 and with that 
for BC = + 0.1195. The aberration is 6 times greater 
for horizontally entering light and 18 times greater 
for light entering the cell under an entrance angle of 


the wrong sign. 


Condition for aberration-free imagery with a 
variable entrance angle. If at any given time we 
choose an entrance angle satisfying condition (81), 
then we know that the function in figure 3 nowhere 
exceeds the value of 0.01421, and if we put down the 
condition 
0.01421 
360B4 


<< O02, 


we get a critical time : 


(84)  B* > 0.001974, B® > 0.04443, B > 0.2108 


after which the aberrations have certainly fallen 
below 1/50 of a light wave in the entire interfero- 


gram. = 
The condition (84) is thus a necessary condition for 
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an aberration-free imagery, and it would also be 
sufficient if the entrance angle were made to vary 
with time according to condition (81). It is of course 
not impossible to do this, but it is a little inconve- 
nient to make such an adjustment before every 
exposure. We will therefore investigate the possi- 
bility of using a constant entrance angle throughout 
the diffusion experiment. It is then quite natural to 
choose that angle which satisfies condition (81) at the 
critical time defined by (84). 

General condition for aberration-free imagery. Inser- 
tion of the value of B from (84) into condition (81) 
gives a numerical value of the entrance angle parame- 
ter: 


(85) C — 0.5669. 


We know that this entrance angle gives an aberra- 
tion-free interferogram at the time defined by (84). 
We now have to prove that the aberration does not 
increase beyond 4/50 in the course of time when C is 
kept at the constant value (85). There is apparently 
a risk for such an increase, since the aberration, 
regarded as a function of time only, was found to pass 
first a minimum, then a maximum, before it falls off 
to zero. The time period between this minimum and 
this maximum is evidently a period of increasing 
aberration and has to be investigated in more detail. 

The time interval for an increasing aberration is 
given by equations (79), in which we now introduce 
the numerical value of C: 


(86) {.0217 f < B < 1.6243 f. 
For f = 0.1298, y = + 1.499, we have the interval 
(87) 0.1326 < B < 0.2108 


from which we can conclude that the aberrations 
will never exceed 1/50 of a light wave outside 
y = + 1.499 after the time defined by (84), if the 
entrance angle is kept constant at the value (85). It 
remains to be seen what happens inside y = + 1.499. 

We then go to the second equation (80) in order to 
investigate the magnitude of the aberration maximum 
that occurs after the tentative critical time (84). After 
introduction of the numerical value of C, we can put 
down a requirement that this maximum be smaller 
than 1/50 of a light wave 


(88) yf < 17.24. 


The function y// increases monotonically from 0 to 
infinity when y passes from 0 to infinity. It acquires 
the value of 17.24 for y = + 1.679, hence the condi- 
tion (88) is satisfied in the region inside y = + 1.679. 
This region overlaps a little the regions outside 
y = + 1.499, which is to be regarded as a good luck, 
since we now can draw the following conclusion. At 
a constant entrance angle, defined by (85), the aber- 
rations are bigger than 1/50 of a light wave at many 
points of the interferogram before the elapse of a 
diffusion time defined by equation (84). During this 
time, each point outside y = + 1.499 has a period 
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of increasing aberrations, with maxima which may be 
greater than 2/50. After the critical time of diffusion 
given by (84), the aberrations lie below 4/50 in all 
parts of the interferogram. It is true that all points 
within y = + 1.499 have a short period of increasing 
aberrations, but these aberration maxima stay below 
4/50. These relationships are also clearly illustrated 
in figure 2. Thus it has been proved that it is feasible 
to use the constant entrance angle given by (85) and 
that it is quite unnecessary to arrange for a varia- 
tion in time according to equation (81), the optimum 
relation for a minimum aberration. 


The Optimum Entrance Angle and the Critical Time 
of Diffusion. If we now go back to the original variables 
8), 0, and ¢t with the aid of equations (65), (71), 
and (72), we get the following expression for the 
optimum entrance angle *: 

we nla 
x De gece ee 
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(89) 


and the critical time of diffusion, at which all aberra- 
tions have decreased to 1/50 of a ight wave, is given 


by the equation 
rna?? 


(90) 2Di = 5? = 0.04443 


With reasonable figures, m = 4/3, 1 = 5461 A, a = 
5 om, and N = 25, we get 5, = 0.0048. This means 
that the light source should be situated about 5 mm 
below the optical axis in the focal plane of a lens 
with a focai distance of 1 meter. The result is quite 
reasonable. 

In view of the fact that the aberration terms con- 
tain the cell thickness in up to the fifth power [ef. 
equation (37)], it is surprising indeed to find that the 
optimum entrance angle and the critical breadth of 
the diffusion boundary (2b) are inversely and directly 
proportional, respectively, to the square root of the 
cell thickness only, whereas the critical time becomes 
directly proportional to the cell thickness in the 
first power. Insertion of the reasonable figures sug- 
gested above also reveals that the critical time is 
surprisingly small. For cane sugar, with D = 5.231. 10-8, 
one gets the figure 109 seconds. For water solu- 
tions investigated by the green light from a mercury 
lamp, equation (90) can be written : 


(91) 0.5280 Dt 


if the total number of fringes is 25, and if t is measur- 
ed in minutes and D in units of the sixth decimal 
place. This leads to the simple diagram shown in 
figure 5. 

It is also of interest to get an idea of the fringe 
density of the interferogram at the critical time. This 
is given by the equation 
On NT fVim VN 

n b 0.2108 \/xa 


* The sign of 39 is opposite that of n/(zx). 


a= 


(92) 
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I'ic. 5. — Graph showing the critical time of diffu- 


sion for different cell thicknesses and_ different 
diffusion coefficients. The graph is calculated for 
water solutions, for the wave-length of the green 
mercury line, and for a concentration of solute 
giving interferograms with 25 tringes. The critical 
time rises with the power 3/2 of the total number of 
fringes. 


At the centre of the gradient, we have f = 0.3989, 
and if we give to the parameters the same values as 
earlier, we get a fringe density of 296 per cm at unit 
magnification, or a fringe separation of 34 ». Thus 
one probably has to use a magnification factor of 
2-3 and fine-grained plates in order to reach the 
resolution of 1/50 of a light wave, which was chosen 
as the basis of our numerical equations. 


Vertically extended light source. It was found in 
the second order theory that a vertically extended 
light source must be expected to give rise to fringe 
blurring unless the Gaussian image plane of the plate 
coincides with the middle of the cell (r = 1/2). This 
conclusion was experimentally verified in reference 2. 

Blurring of fringes arises from aberration terms 
containing 3) as a factor. This is natural, because 
light from the upper and lower ends of an extended 
light source must necessarily strike the entrance wall 
of the cell under different entrance angles. If 38, 
is retained to denote the entrance angle of light 
coming from the centre of a vertically extended light 
source, then every point of the cell will receive light 
pencils with all entrance angles lying between two 
limits, say 8) + ¢, where 2c is the angular extension of 
the light source when viewed from the cell. It is then 
easily realized that aberration terms containing 8) 
in the first power will, in addition to the aberration 
they cause for a point source, give rise to a blurring 
which is symmetrical towards both sides of the fringe, 
that is, to a blurring not accompanied by any fringe 
displacement. On the other hand, terms containing 
the squared entrance angle must give rise to an asym- 
metric blurring and to a displacement of the centre 
of the fringe. 
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The most important blurring aberration is of course 
that already derived in the second order theory. This 
term is not present in equation (60) because of the 
direct transformation to the coordinate #,, and 


“m) 
because the equation was derived on the assumption 
of a point-shaped light source. In reference 1, however, 
the following equation was derived for the < value cor- 
responding to complete blurring : 

(93) e a? n’ (1—2r) = ma 

Automatic guard against aberrations in interfero- 
grams. Half-blurred fringes are apparently obtained 
for the half of this ¢ value. After introduction of the 
variables y and B, we thus get the following condi- 
tion for half-blurred fringes : 


8B xm 
af v/a NA 


The minimum permissible value of B is 0.2108, and 
the value of f at the centre of the boundary is 0.3989. 
tlence the fringes at the centre of the interferogram 
will be half-blurred at the critical time if the angular 
extension of the light source satisfies the equation 


(94) << 


(95) = 


Apparently it is a clever arrangement to use such 
an e¢ value since the blurring of the central fringes 
will then automatically prevent a meaningless 
evaluation of interferograms exposed before the cri- 
tical time of diffusion has elapsed. In other words, 
by using a light source with the vertical extension 
given by (95), the fringes will grow sharp and well- 
defined at the same time as the aberrations are fading 
out to insignificance. 

With the same numerical figures as used earlier 
in this article in order to exemplify applications, we 
get 2 ¢ = 0.0047, which means that a 4.7 mm long 
vertical slit should be used in the focal plane of a 
lens with a focal distance of 1 meter. 

The additional blurring due to third-degree terms 
containing the entrance angle as a factor need not be 
considered. 


2. The absolute aberration function in uncollimat- 
ed light. 


Additional aberration terms in the case of uncolli- 
mated light. Since the entrance angle 3, appears only 
in third-degree terms, it is permissible to use the 
value: 


t—4q x 
8, = + ~— + p 
p P 


(96) 


instead of the exact value given by equation (5), 
because the quantities (x,,—w,)/p and ra/p are of the 
second and first degrees, respectively. In equation 
(96) p will then be regarded as a divergence para- 
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meter, while e is the off-axis parameter. The latter is 
identical with the entrance angle at the centre of the 
diffusion boundary and thus takes over the role hitherto 
played by 8). 

On introduction of (96) into the absolute aberra- 
tion function, given by the last four terms of equa- 
tion (61), the terms containing 9 but not p become iden- 
tical with the function (63) already treated, with 0 
instead of 8). Then there remain four terms contain- 


ing p: 


Qn" 2 8p" 
ans atxn" ep 
97 Lp 6) —— = aN 
2H) oh) et Pe) 24m2p?r 12m?pr 
aten'n" ain’? 
24m? pr 72m*pr 


After introduction of the variables y, B, and C with 
the aid of equations (64), (66), (71), and (72) (with ¢ 
instead of 8), and after elimination of p in favour of 
the new dimension-less parameter P defined by the 


equation: 
1 Pm 
98 maps 
( ) p a\ N 


the function (97) assumes the simpler form : 
Pey?f  PCy?f 
24 2B 
Ro ja ei) 
ihe £127. 


(99) @ (y,B,C,P,) = 


or with the numerical values of B and C according 
to equations (84) and (85) : 


(100) Q(P,y) = — 0.0417 P2y%f + 0.2241 Py?f + 
+ 0.3126 Pf? (1—3y?) 


The permissible divergence or convergence of a 
light beam. The course of the divergence aberra- 
tion (100) through a diffusion boundary is given in 
figure 6 for P = 0.4, which is the numerical value of 
the divergence parameter that gives rise to aberra- 
tions touching our tolerable limit, 1/50 of a light 


RESOLVING POWER OF INTERFEROMETRIC METHOD 


Fic. 6. The contribution to the absolute aber- 
ration A which is due to non-parallel light, for a 
divergence parameter P = 0.4. The figure shows this 
divergence to be the maximum permissible. 
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wave. The same negative value, which corresponds to 
convergent light through the cell, is also useful. As 
a matter of fact, P = — 0.4 gives a curve centrosym- 
metric with respect to that in the figure, since a 
sign reversal of both P and y in (100) changes the 
sign of each term. 

We can thus formulate the following condition for 
the divergence parameter in the case of uncollimated 
light : > 

(101) pele ee 

mp 


which can also be written in the form 


(102) 


With the same numerical figures as applied earlier, 
one finds that a/p should be smaller than 0.11, which 
is also a reasonable figure. 


3. The relative aberration function in collimated 
light. 
This function is given by the equation 


2 ° 

by q 85 adn’ ata"? 
(103) R(z,t, 89) = —-~ + ~—— + —. 
2m? 2m? 6m* 


After introduction of the variables y, B, and C, it 
takes the form 


(104) 


/n “1 es 
Rly,B,C) = WA jo: ah gad 


2ma Bo 23Ree 
With a fixed value of C, this relative aberration can 
be regarded as a function of one single variable, f/B. 
This function has a minimum and increases in both 
directions therefrom. Consequently the highest value 
this function can acquire is either at the lower or 
at the upper limit of the interval within which //B can 
vary. These limits are 0 and 0.3989/0.2108 = 1.8923, 
respectively, and the latter limit gives the higher 
value of AR. Hence this aberration satisfies the 
inequality (for C = — 0.5669) : 


0.22112. /N 


mea 


(105) Ry, B, OQ < 
With the very unfavourable figures of V¥ = 100 and 
a = 1 cm, this gives R = 0.9-10-4, which is a relative 
error without importance. It may be concluded that 
the relative aberration never requires consideration 
under conditions making the absolute aberration 
insignificant. 

It should be noted that negative C values (if f 
positive) have the same favourable effect on this 
error as it was shown to have on the absolute aber- 
ration. If C were allowed to change sign, all three 
terms of R would get the same sign. 


4. The relative aberration function in uncollimated 
light. 
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Introduction of (96) into (103) gives the following 
three terms containing the divergence parameter p- 


If 


axn 
2m*p 


. x2 . OL 
HG Oe a 0) We 2m2p? | 2m*p 


(106) 


The first two terms have their origin in the cosine 
error due to obliquely incoming light. This is the only 
aberration that actually increases in the course of 
time. This is evident from the fact that the « interval 
covered by the diffusion boundary grows larger and 
larger indefinitely as the diffusion process progresses. 
Thus the first term in the above equation, especially, 
increases with time rather rapidly and without 
limit. 

An upper bound to this aberration can best be obtain- 
ed by insertion of plausible figures in each term 
separately. There is no reason for making a diffu- 
sion cell longer than 6 cm, thus x is always smaller 
than 3 cm. Likewise, one can say that a plausible low- 
er limit for p is 50 cm (cf. equation 102). With m = 
4/3, one then finds that the first term of (106) stays 
below 10-8, which is at the border of a tolerable 
error in precise measurements. 

In the second term, we will assume pe to be smaller 
than 0.01, and make the same assumptions regard- 
ing x and p as above. We then find that the nume- 
rical value of this term never exceeds 1.7-10-4. 

The third term actually contains the second deri- 
vative, since it can be written in the form 


(107) aan’ abn’ —a(An)f' » a (An) yf 
2m2p  —-2m?p 2m?p 2m2p — 
The maximum value of f' = — yf is 1)\/2ne = 0.242. 


With the same numerical values as earlier for m 
and p, we thus find this term always to be smaller 
than 4-10-5, if An is smaller than 0.01. 


VII. Discussion. — Mathematical considerations. — 
Although only elementary mathematics is involved 
in this treatment, the problem must be said to be 
unusually complex. The analysis is connected with 
an extensive algebraic labour. The risks of making 
algebraic mistakes are numerous, and the computa- 
tions have to be carried out with a scrupulous care. 
It should be sufficient for the interested reader, how- 
ever, to get a clear understanding of the principles 
of the mathematical methods. The correctness of the 
equations is granted by the fact that the same results 
have been gained by independent mathematical 
methods. It would lead too far to describe in detail 
other possible ways of computing the third order 
aberrations, but one hint will be given. The second 
order approximations for the exit coordinate and 
the exit angle can be derived by the method of 
undetermined coefficients. It one puts, tentatively, 


(108) 2,==a,+ Aa, +B (An) + Co, +Da, (An) + 
+ E (An)? 
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it follows because of equation (13) that the exit 
angle is obtained by differentiation with respect 
toa: 


dA dB OG 5 
109 = c= An) 4 s 
ee) ae Paar Pe a aaa 
dD dE 
da °° Ces da eel 


If these expressions are introduced into equation (7), 
and if n(x,) and the cosine functions are developed 


into their second order approximations, an identity 
is obtained which in turn gives rise to a series of 
differential equations. These can easily be solved for 
the coefficients A, B, C, D, and EF. The result thus 
obtained is identical with equations (21) and (26). 

It is also possible to check the correctness of the 
equations by insertion of the different variables into 
the exact equations (6) — (10). The third degree approx- 
imations of these equations then must reduce to 
identities 0 = 0. One such check will be demonstrat- 
ed in Appendix 1. 

The author has attacked the theory of the light 
path through a diffusion boundary on two earlier 
occasions [5, 6]. The first trial was made in 1940 
using very simple mathematical means, while the 
later one was more elaborate. Both these works 
pertained to the Schlieren method, so the results 
of the present theory have to be differentiated with 
respect to x before a comparison can be made. They 
are open to criticism from both mathematical and 
optical points of view, but, nevertheless, it is interest- 
ing to observe that many of the results now reached 
on firm and sound mathematical and optical bases 
were actually anticipated already in the 1940 article. 
Thus it was concluded already there that the WIENER 
skewness aberration could be removed, although 
a wrong numerical value of the cell defocusing was 
reported. The flattening of the top of the gradient 
curve, the correct value of which is obtained by 
differentiation of equation (63) with respect to # and 
by putting a. 10: (77 =.0): 


dA(0,t, 8) — a885n""(0) 
(110) ae 
daz 24m? 
atdon'(O)n""(0)  an'(0)n'""(0) 
24m? 90m? 


was reported in the 1940 article to have the figure 72 
instead of 90 in the denominator of the last term 
(only the case of horizontally entering light, 8 = 0, 
was investigated). This figure was the result of a 
guess between two limits, 48 and 6, given by the 
mathematical analysis. In the paper of 1946, the cor- 
rect value of this term was found (p. 79, equation 79). 
As we have seen in this article, the Wiener skew- 
ness and the aberration (110) are the two most im- 
portant aberrations. A critical time condition was 
also given in reference 6 (p. 80, equation 103), accor- 
ding to which this time appears to be proportional 
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to the square of the cell thickness. This discrepancy 
relative to the result reported here is, however, only 
apparent, since our parameter N includes the cell 
thickness in the first power. 

The fact that a defocused cell is most favourable 
even in the cases of finite entrance angles and uncol- 
limated light was not revealed in these earlier works. 

The fact that different aberration terms of the third 
degree have been found to counterbalance each other 
to a great extent is very valuable, since it makes it 
possible to bring down the aberrations quite appre- 
ciably. However, it also involves a risk that fourth- 
degree terms, if they do not balance each other in 
the same favourable way, are not appreciably smaller 
than the algebraic sum of the third-degree terms. 
The theory retains its essential validity, however, 
even if the fourth-degree terms rise to the same order 
of magnitude as the minimized sum of the third- 
degree terms. The results gained cannot be quest- 
ioned unless one assumes the fourth-degree terms 
appreciably to surpass the third-degree terms in 
magnitude. Such an assumption appears to the author 
as highly improbable, if not impossible. Anyhow, 
experimental proofs of the validity of these results 
must be regarded as highly desirable. 


Optical considerations. — The author has been 
criticized regarding the work in reference 6 on this 
topic for the reason that the Schlieren method can- 
not on the whole be treated by ray-optical means. 
This criticism was justified, and it was also shown 
through the excellent works of KEGELES and Gos- 
TING [7], LoneswortH [8], CouLson et al. [9], 
ADLER and Briancuarp [10], and Gostina and 
OnsaGcer [11] published shortly after reference 
6 that the various modifications of the Schlieren 
method, the scale method, and the Gouy interfe- 
rence method are afflicted with an aberration at the 
top of the gradient curve that could only be account- 
ed for by a wave-optical theory. The objects of the 
present theory are interferometric methods charac- 
terized by optical imagery of the cell, and the diffi- 
culty discussed above does not enter the picture. It 
should be pointed out in this connection that inter- 
ferometric methods for recording the gradient curve 
are available [12, 13, 14, 15], and that these methods 
are also free from the wave-optical aberration dis- 
cussed above. 

The analysis presented here has been guided by 
experimental-optical considerations, whereby care 
has been taken to avoid unnecessary assump- 
tions based on conventional habits and precon- 
ceived opinions regarding optical adjustments. This 
principle has resulted in the introduction of a number 
of parameters, the numerical values of which are 
open to free choice, each parameter having an expe- 
rimental-optical counter-part. These parameters are 
the cell position parameter, the entrance angle para- 
meter, the divergence parameter, the parameter 
related to the vertical extension of the light source, 
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and the parameter pertaining to the plane in which 
the object is placed in the calibration of the plate 
coordinate in terms of the cell coordinate. The last 
parameter was never introduced, because its opti- 
mum value was already known from the second 
order theory. That was the case with the cell position 
parameter too, but this parameter was introduced 
here again in order to show that its optimum value, 
r = 2/3, also minimizes the third-degree aberrations. 

The introduction of so many parameters of course 
contribute to make the handling of the mathematics 
extremely complex, but only in this way has it been 
possible to reach a number of valuable conclusions. 
Only in this way could it be revealed that a certain 
defocusing of the middle of the cell has a very pro- 
nounced and favourable effect in that it completely 
removes the important aberration corresponding 
to the well-known Wiener skewness of the gradient 
curve. Likewise, the far from self-evident recommen- 
dation of measuring the photographic enlargement 
from the middle of the cell, in spite of the defocusing 
mentioned, could not have been found by a less 
elaborate procedure. 

The entrance angle has hitherto been regarded as 
more or less arbitrary, and horizontal, collimated 
light is in general used. It has now been shown, 
however, that a certain entrance angle of a sign 
opposite to that of the refractive index gradient gives 
much smaller aberrations than horizontal light, and 
appreciably reduces the time required for the aber- 
rations to fade out to insignificance. The equation 
given for this optimum entrance angle is easily appli- 
cable to experimental work since it only contains 
known apparatus constants and data available from 
the preparation of the experiment (the total number 
of fringes — proportionnal to the concentration and to 
the cell thickness). 

The dimensions of the light source have also been 
regarded as more or less arbitrary hitherto, provided 
that they satisfy the requirements of the interfero- 
metric optical system without regard to the diffu- 
sion cell. As is well known, some such systems require 
a point source, others can be operated with a slit- 
shaped source, while others again may be run with 
a circular or rectangular source. In the latter two 
cases, the stop at the light source has been used 
merely in order to regulate the light intensity to a 
proper value. In the previously published second- 
order theory, it was found that a vertically extended 
light source must give rise to fringe blurring as soon 
as the middle of the cell is defocused. Thus the defo- 
cusing necessary for abolishing the Wiener aber- 
ration has to be accompanied by a suitable reduction 
in the vertical extension of the light source. It has now 
been shown in the third order theory that this fringe 
blurring actually can be used to advantage. By a 
proper adjustment of the vertical extension of the 
light source, it is possible purposely to produce blur- 
ring of the central fringes during an initial time period 
while the interferograms are not yet aberration-free. 
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The necessary light intensity can be obtained by a 
proper choice of lamp and by selecting an interfero- 
metric set-up that allows the use of a horizontally 
extended light source. 

The realization of such an automatic guard against 
a meaningless evaluation of interferograms not yet 
free from aberrations necessitates exact equations for 
the fringe blurring and for the time required for 
aberration-free imagery. Both these equations are 
given. They are also combined into an explicit and 
easily applicable equation for the optimum vertical 
extension of the light source. 

The critical time period, ¢. e. the time required for 
the aberrations to decrease below the resolving power 
of the optical system, has been found to be surpris- 
ingly short if the tricks recommended in the optical 
adjustment are applied. The diagram in figure 5 gives 
an idea of the order of magnitude of this critical time 
for different diffusion coefficients and cell thicknesses. 

It has also been shown that precise measurements 
can very well be carried out using uncollimated 
light. There is no optimum value of the divergence 
or convergence parameter, because the additional 
aberrations caused by uncollimated light have the 
same sign through the whole diffusion boundary, 
while the aberrations in collimated light have oppo- 
site signs on both sides of the starting boundary, 
An equation for the maximum allowable divergence 
or convergence of the light is given. 


Physico-chemical considerations. — Experimental 
workers in diffusion are familiar with two different 
deviations of refractivity gradient curves from the 
shape of the Gaussian error function, which is obtain- 
ed in the case of ideal diffusion. One deviation is a 
skewness, and its physico-chemical background is a 
concentration-dependent diffusion coefficient. The 
other deviation consists in an abnormal height of 
the curve. This deviation is the result of polydisper- 
sity or of impure preparations. The study of these 
deviations from the normal error curve is of great 
importance, since they involve possibilities of measur- 
ing concentration dependence and_ polydispersity 
if and when reliable theories of these phenomena 
become available. 

It is interesting to note that the second and third 
order aberration theories have revealed the presence 
of two optical aberrations of the same nature as those 
mentioned above with physico-chemical origin. 
It is quite natural, then, that it is impossible to make 
correct physico-chemical interpretations of non-ideal 
diffusion curves unless they are optically aberration- 
free, or corrections can be applied. 

The possibility of removing the optical skewness 
aberration completely is consequently a matter of 
great physico-chemical interest. Thus it has been 
found recently by THompson [16] that the skewness 
in early exposures differs from that in later expo- 
sures for material with a concentration-dependent 
diffusion coefficient, a variation that could be explain- 
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ed by taking the Wrener skewness into account. 

The absolute aberration function (63) results in a 
flattening of the gradient curve, that is, to an effect 
opposite to that of polydispersity. This depression 
of the n'(z) maximum is given by equation (110). 
An interferogram photographed before the aberra- 
tion (63) has decreased enough may therefore give 
rise to an erroneous interpretation of the polydisper- 
sity of the material. The direction of this misinter- 
pretation is such as to indicate too small a hetero- 
geneity. 

In case of a pronounced concentration dependence, 
it is customary to measure the diffusion coefficient 
at a series of different absolute concentrations, using 
a small and constant concentration difference bet- 
ween top and bottom solutions in the diffusion cell. 
In order to retain a high accuracy in such measure - 
ments, it is desirable to increase the cell thickness 
in order to keep the total number of fringes at a 
suitable value (25-50). In this connection it is of 
great interest to note from the present theory that the 
critical time period rises with the first power of the cell 
thickness only, although the aberration terms them- 
selves contain the cell thickness in the fifth and lower 
powers. This fact, together with the favourable indi- 
cations regarding the critical time period shown in 
figure 5, promises very good conditions for measuring 
diffusion coefficients at extremely low concentration 
differences. Its lower limit depends probably more on 
the temperature regulation device than on the accu- 
racy of the optical recording system. The lower the 
concentration difference, the greater the sensitivity 
of the boundary to liquid circulation below and on 
top of it. 


Appendix I. — Here we will demonstrate in detail a 
cross-control of the correctness of two equations of 
the theory. We choose the equations (25) and (29), 
and the check will be made using equation (10). 

It is necessary during such a computation always 
to be clear about the degrees of the different terms. 
Otherwise much work in vain is carried out in car- 
rying along terms of the fourth and higher degrees, 
with which we are not concerned in this theory. Thus, 
the square of the external exit angle is found to be : 
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On the left-hand side of equation (10), it is. suffi- 
cient to use the first order approximations of the sine 
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functions, that is, to put the sines equal to their 
angles, because the third degree terms on squaring 
give rise to terms of at least the fourth degree. Thus 
we get, with the aid of (111), 
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This is the third order approximation of the left- 
hand side of equation (10). 

If we now turn to the right-hand side and use equa- 
tion (25), it is again advantageous to distinguish 
sharply between first-and second-degree terms. If we 
denote them by subscripts 1 and 2, respectively, we 


can write 
(113) n(%,) = n (2,) (@,—2%,), 7° (a). + 
l " 
(Z,—2p)a N (24) ) i TT] (x) 


Here the first term on the right hand side is of the 0: 
th degree, the second is of the second degree, and the 
third and fourth terms are of the third degree. Hence 
first-degree terms are not present. When we now have 
to square this equation, and want to retain only 
terms of up to the third degree, there only remain 
the square of the first term and the double products 
of this term and each of the other three terms : 
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In this equation we put mo = m + (no—m), develop 
the square, and discard the terms containing the 
square of (7 —m) since they are of the fourth degree. 
We then obtain 


(115) n* (%,) =n? (a) = 2mas,ni + an’? + 
uu 4 ora bath, on, 
+ maro.n, + | 
3 3m 


Comparison with equation (112) shows that equa- 
tion (10) is in complete agreement with the second 
order approximations (25) and (29). 


Appendix Il. The aberrations caused by plano- 
parallel plates. — In figure 7 there is a plano-parallel 
plate of thickness t and refractive index m, situated 
between the planes z = —t(1—1/m) and z =t/m. The 
plane through the x axis perpendicular to the plane of 
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the paper is then to be regarded as the Gaussian 
image plane of the left-hand plane surface of the 
plate. We will regard a light pencil entering the plate 
under the entrance angle 8 on the z axis. Within the 
plate, it makes the angle e with the axis. At B the 
pencil leaves the plate, and when traced rectilinearly 
backwards from the exit point we arrive at the point 


lic. 7. — Optic effect of a plano-paraltlel plate. The 
v axis is the Gaussian image plane of the left-hand 
surface of the plate of thickness ¢. The origin, O, is 
thus the optical image of the point A, and the 
aberration OA/ is calculated in the text, as well as the 
optical path length between A and A’. 


A’ in the Gaussian image plane of A. In the first 
order approximation, the distance OA‘ = Ax is equal 
to zero, but we will now calculate the third order 
approximation as a function of 8, m, and ¢t. We will 
also compute the optical path length between the 
object A and its image A’ as a function of the same 
variables. This will give us a quantitative measure of 
the deviations from true Gaussian imagery through 
plano-parallel plates. 


The vertical shift of the Gaussian image. The magni- 
tude of Ax can be taken directly from the figure 


(116) Av = ¢t tan 9 — cae tan § 


m 


and the third-order approximation of SNELL’s law of 
refraction is 


38 m—1 
OSS 3, 
m 6m3 


(117) 


After the tangent functions in (116) have also been 
developed into their third order approximations, 
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equation (117) is introduced into it, which gives the 
following equation : 


m?—1 


118 
(118) oe 


A (1 

The step from equation (59) to (60) in the present 
analysis consisted in the removal of a planoparallel 
plate of thickness p and refractive index m on the 
assumption that the imagery through this plate were 
strictly Gaussian. Equation (118) now gives us the 
possibility of investigating the order of magnitude 
of the deviation from Gaussian imagery through 
such a plate. Since p varies between, say, 50 cm and 
infinity, this deviation appears, at the first sight, to 
increase indefinitely on passing over from uncollimat- 
ed to collimated light. This is incorrect, however, 
because the two coherent pencils have nearly the same 
direction. The error that comes into play is thus only 
the difference in (118) for the two angles 8, and 8,, and 
we get with the aid of equation (47) (7 is put = 2/3): 


eo es 


(119) 


4m* 


With the aid of relations given earlier, this can be 
brought into the form 


(m?—1) C? FN 3? 


4Bat?ml 


(120) NY = 


and the error has its highest possible value at the 
critical time, B = 0.2108, with C = — 0.5669, and 
svrehepe = 03989. With NV =925) a= 5 em, m= 4/3, 
and 2 = 5461 A, the numerical value of this error 
is found to be as low as 11.6 Angstrém units. 

On the camera side of the cell, we have to take 
into account the second cell wall and, possibly, a 
thermostat window and a certain thickness of ther- 
mostat water. The optical imagery through these 
plano-parallel plates has also been assumed to be Gauss- 
ian. We can assume the added thickness to be less 
than 10 em of water, and the two angles we have to 
deal with are 8, for the one and 8 + an’ for the 
other pencil. One easily finds that this error is of the 
same order of magnitude as that on the light source 
side of the cell. 
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The error in optical path length. The optical path 
length from A to A’ in figure 7 is given by the equa- 
tion 


(121) mt t 


Ss = 


COS © m cos 8 


The cosine functions here have to be approximated 
to the fourth degree: 


1/eos 7 = 1 + 2? /2+ 5 x4 /24. 
We then obtain, with the aid of (118), 
Panes / 294 
(122) es t (m?—1)_ pee 38 7 
m 8m? / 


Here again, we are only concerned with the difference 
in path length for the two pencils. Hence we have, on 
the light source side of the cell, with t = p, and the 
angles equal to 8) and 8, : 


5 3 

(m?—1) 8,a?n' 
AS 
4m* 


Comparison with equation (119) reveals that this error 
is still one order of magnitude smaller, consequently 
far below the measurable range. 
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The determination of the absorption coefficient of luminescent microcrystals 


Gy. GERGELY 
Research Institute for Telecommunication, Budapest, Hungary. 


Summary. — The absorption and luminescence emission of single Zn-type microcrystals was investigated. A photoelectric microscope 
photometer was built which consists of a ZE1ss microscope and a 931 A multiplier phototube. Evaluating the results of absorption 
measurements, neither a plane-parallel, nor a spherical crystal model furnished reliable results for the absorption coefficient. The 
application of an empirical model, based on the non-absorbed radiation (5 770 A) made possible the determination of the ultra- 
violet absorption coefficient by comparing the transmission of single microcrystals for different wavelengths. The method was 
experimentally checked on K,Cr,0, crystals. The ultraviolet absorption coefficient of luminescent microcrystals may also be 
determined from the dependence of the light emission on the particle size. 


Sommarre. — On éludie ici Vabsorption et ’émission luminescente d’un microcristal simple de Zn-type. Un microphotométre photo- 
électrique a été construit dans ce but, & partir d’un microscope Zetss et d’une cellule photomultiplicatrice 931 A. L’interprétation 
des mesures d’absorption n’a donné de résultats valables ni pour un cristal plan paralléle ni pour un cristal sphérique. L’utilisa- 
tion d’un modele empirique, basé sur la radiation non absorbée (5 770 A) a rendu possible la détermination de coefficient d’absorp- 
tion de Vultraviolet, par comparaison avec la transmission de quelques microcristaux pour différentes longueurs d’onde. La 
méthode a été appliquée expérimentalement pour des cristaux de K,Cr,0,. Le coefficient d’absorption de Vultraviolet de micro- 
cristaux luminescents peut aussi étre déterminé en relation avec l’émission lumineuse de la particule. 


ZUSAMMENFASSUNG. — Die Absorption und Lumineszenzstrahlung eines einzelnen Mikrokristalls vom Zn-Typ wurde untersucht. 
Hierzu wurde ein elektrisches Mikrophotometer aus einem ZE1ss-Mikroskop und einem Multiplier 931 A zusammengestellt. 
Bei der Auswertung der Absorptionsmessungen ergeben sich aber weder fiir eine planparallele Platte noch fiir einen sphdrischen 
Modellkristall zuverldssige Werte des Absorptionskoeffizienten. Erst die Verwendung eines empirischen Modelles, das von der 
nicht absorbierten Strahlung (5 770 A) ausgeht, machte es mdglich, den Absorptionskoeffizienten im Ultraviolett zu bestimmen, 
indem man die Durchldssigkeit einzelner Mikrokristalle fiir verschiedene Wellenldngen vergleicht. Die Methode wurde experimen- 
tell an K,Cr,O,-Kristallen erprobt. Der Absorptionskoeffizient im Ultraviolett kann bei einem lumineszierenden Mikrokristall 
auch aus der Abhdngigkeit der Lichtemission von der Teilchengrésse ermittelt werden. 


A considerable proportion of luminescent mate- 
rials are microcrystals of 1-30-100 micron particle size. 
In view of this, the absorption coefficient cannot be 
determined by the simple method as in the case of 
macroscopic crystals. Although the problem is a clas- 
sical one, its solution has only recently been achieved. 
The coefficient of absorption is sometimes determined 
by statistical methods, from the coefficient of diffuse 
reflection, as it was described by Bopo [1], Joun- 
son [2] and Antonov-RomMAnovsky [3], or by the 
microscopic method, developed by BrumBErG and 
PEKERMAN [4]. 

We do not wish to give a detailed discussion of the 
different methods, but some difficulties of the statis- 
tical method must be mentioned. Everybody using 
statistical methods has employed models essentially 
simpler than reality (plane-parallel crystals). The 
absorption coefficient has been expressed in involved 
formulae, and the formulae — excepted Bono’s — 
had an unknown parameter, i. e., the refractive index. 
The formulae also had a term, ul, product of the 
absorption coefficient » and the light path J. The latter 
corresponds to the particle size only in exceptional 
cases. For the most part, the light rays are repeatedly 
totally reflected within the crystal [3]. For this rea- 
son, statistical methods give a smaller value for u than 
the real one. 

Besides the theoretical difficulties, there is a prac- 
tical difficulty, too : the measurements of the diffuse 
reflectance must be carried out on samples having 
homogeneous particle size. In practice, the homoge- 


neous fractions are produced by sedimentation [1 | 
and a considerable quantity of material (50-100 g) is 
necessary for this purpose. A further difficulty lies in 
the fact that repeated sedimentations cause surface 
deterioration, and darken the ZnS type crystals. The 
darkening strongly affects the diffuse reflection. 

The advantage of the statistical methods lies in the 
simplicity of the reflection measurements. Although 
the microscopic method is experimentally more diffi- 
cult, its use is justified by its advantages. The micros- 
copic method was first employed by BrumsBere and 
PEKERMAN [4]. They used a quartz microscope spec- 
trograph for the determination of absorption spectra 
of microcrystals. They photographed the light trans- 
mitted by selected plane-parallel crystal samples. The 
original BRUMBERG-PEKERMAN method offers some 
difficulties. It can only be applied for strictly plane- 
parallel crystals, whereas, the ZnS type microcrystals 
have approximately spherical or ellipsoidal shape. 
BruMBERG and PeKeRMAN measured the absorption 
spectra only up to 4000-4500 A. Otherwise, the pho- 
tographic procedure is not possible. 

To eliminate the difficulties of the original Brum- 
BERG-PEKERMAN method, a photoelectric microscope 
photometer and a new evaluation of the experimental 
results was developed. 


Experimental arrangement. — Figure 1 shows the 
experimental arrangement. The radiation of a high 
pressure 125 W mercury lamp is projected on the 
microcrystal to be examined by means of a mirror and 


vol. 8, n° 4, Dec. 1956] 


ABSORPTION COEFFICIENT OF LUMINESCENT MICROCRYSTALS 


185 


Fic. 1. — Experimental arrangement. Notation on the figure : 1, mercury lamp ; 2, mirror ; 3, condenser lens ; 
4, diaphragms: ; 5, microscope stage ; 6, slide ; 7, microcrystal; 8, objective; 9, binocular attachment ; 10, 


eyepiece ; 11, eyepiece micrometer ; 12, diaphragms 


; 13, prism ; 14, light modulating sector ; 15, motor; 16, 


filter; 17,/multiplier phototube ; 18, load resistance ; 19, electronic voltmeter ; 20, high voltage supply. 


a condenser. The light rays are limited by suitable 
diaphragms. The microcrystals are sedimented on a 
glass slide and immersed in vaseline oil. 

Light leaving the crystal is collected by the objec- 
tive. The microscope has a binocular attachment, divi- 
ding the light into two separate beams ; the image of 
the crystal is observed by the left eyepiece provided 
with an ocular micrometer, in order to determine the 
particle size and to adjust the object on the optical 
axis. 

A real image of the object is produced on the cathode 
of a multiplier phototube by the right-hand beam, 
diaphragms being provided to eliminate stray light. 
The upper diaphragm selects a 4.5 micron diameter 
part of the microcrystal on the image, so that light 
from a small area only of a larger (10-40 micron) crys- 
tal reaches the photocathode. 

The multiplier phototube was a selected sample of 
type 931 A. The light entering the photocell was modu- 
lated by a rotating sector. The photocurrent produced 
an alternating voltage on the load resistance, measu- 
red by an electronic voltmeter. The photocell was 
combined with suitable optical filters (SHort interfe- 
rence, Scuotr Glass, Jena). 

The total magnification of the optical system was 


160. 


Since a quartz microscope was not available, mea- 
surements were carried out only up to 3 650 A. 

The measurements were performed as follows : the 
individual microcrystals were introduced onto the 
optical axis. Using different optical filters, the radia- 
tion escaping both from the single crystals and from 
the ‘‘ background ” was measured. The intensity of the 
‘‘ background ” radiation is equal to the radiation 
reaching the crystals. The measurements were carried 
out with many single crystals of different sizes. 


The evaluation of the experimental results. — 
The photoelectric microscope photometer furnishes 
the following data : the particle size of the individual 
crystals, the intensity of radiation leaving the crystal 
and the intensity of the ‘‘ background ” radiation. 
The 3 650, 4 046-78, 4 358, 5 460 and 5 770-90 A mer- 
cury lines were sufficiently intense to be detected des- 
pite of the very strong attenuation in the optical sys- 
tem. 

To evaluate the experimental data, three different 
models have been considered. 


a) The plane-parallel model. — The plane-parallel 
crystal model was first used by Bopo in his statistical 
method. Although good results were obtained by the 
statistical method, the plane-parallel model cannot be 
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used in case of individual crystals, because their shape 
differs considerably from this form and absorption 
measurements on single crystals an incompatible with 
the plane-parallel model. ZnS microcrystals have 
negligible absorption of the 5 460 and 5 770 A mercury 
lines [7]. In the case of non-absorbed radiation, the 
transmission of a single crystal would be given by : 


1 —o (/n—1? 
— 
Tina’ \n44/ 

where n is the refractive index. For ZnS crystals 
immersed in vaseline oil, n Y 1.5, giving 


Tess O92 


{ss 


independently of the particle size. The experiments 
disproved this supposition. 

b) The spherical model. — ANTONOV-ROMANOVSKY 
used the spherical crystal model in his statistical 
method. Examining ZnS type microcrystals by the 
microscope, their shape was found to be approxima- 
tely spherical or ellipsoidal. When parallel radiation 
reaches a spherical crystal, the particle acts as a sphe- 
rical lens. Mecke [5], Blumer [6] and others dis- 
cussed the optical properties of spherical particles. 
Their problems, in some aspects, differed of those pre- 


Fie. 2. — Ray tracing in a spherical crystal, 
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sented by our work. Using the optical system of 
figure 1, the rays leaving an area of 4.5 micron diame- 
ter of the object reach the photocathode. Figure 2 
shows the ray traces within the spherical crystal. The 
radiation transmitted by a spherical crystal was cal- 
culated by integrating the Fresnet-formulae on the 
surface of the sphere. The calculations are of interest 
for non-absorbed radiation, because they can be 
checked experimentally. It was found that the trans- 
mission of a spherical particle is given by : 


fie LEI 


denoting by / the diameter of the sphere. G(/) is a com- 
plicated integral based on the Fresnet-formulae. 

The results of experiments did not agree with the 
theoretical calculations. The light collecting effect of a 
spherical crystal predicted by calculations was found 
to exist qualitatively, but its magnitude was less than 
the calculated values. No pronounced correlation was 
found between /* and the transmission of the crystals. 
The discrepancies between theory and experiment 
may be explained by small deviations from the ideal 
spherical shape and by the sub-microscopic roughness 
of the crystal surfaces. The calculations are valid only 
for optically polished spherical surfaces. 


c) The empirical model. — The use of non-absorbed 
radiation makes possible the determination of the 
ultraviolet absorption coefficient in the following 
manner. When a crystal was irradiated by non-absorb- 
ed radiation, it either collected the light rays because 
of its spherical lens properties, or else scattered them 
because of its irregular shape. Denoting by J, the 
intensity of the background radiation, and by J the 
intensity escaping from the crystal, the quantity f, 
is the empirically determined shape factor of the par- 


ticle : 
(RESO Meme 


The quantity f, may be supposed to depend on the 
shape, size and surface of the crystal. It may be assu- 
med that it does not vary with the wavelength ; ‘this 
is only an approximation, since f, depends on the 
refractive index and consequently on the wavelength. 
The small size of the particle (10-40 microns) permits 
the approximation. 

In the case of absorbed radiation, the intensity leav- 
ing the crystal and reaching the photocell, is given by 


Ly Slot ea! 


Lis the light path in the crystal. Because of the approxi- 
mate spherical shape, / was taken equal to the diame- 
ter. In the case of particles of ellipsoidal form, | was 
taken equal to the smaller diameter of the crystal. 
Crystals deviating strongly from spherical or ellip- 
soidal shape were excluded from the measurements. 

Table I shows the experimental results for a 90 ZnS- 
10 CdS-Cu phosphor. It may be seen that the absorp- 
tion coefficients determined on different crystals have 
the same magnitudes. Table I contains f, for two wave 
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TABLE I 
l x lo lo iS 
5770-90 A | 5460 A | U (3650) | uw (4060) | w (4358) 

1.325 1.31 1.17 627 137 107 
SS 1.55 1.58 583 334 178 
14.9 1.19 1.28 678 235 114 
14.9 132 1.44 277 124 124 
14.9 0.97 0.94 624 218 168 
15.6 1.33 Ly 640 250 157 
15.6 1.17 1.09 782 286 308 
16.9 0.885 0.885 438 93 77 
16.9 0.72 0.73 418 7 77 
17.6 1.64 1.88 765 762 131 
18.2 0.815 0.790 335 256 49 
20.3 1.56 1.78 582 221 160 
20.3 1.41 1.68 313 191 40 
21.6 Datel 2.41 540 192 88 
21.6 0.483 0.500 304 172 40 
22.3 JES LD 502 238 73 
D316 0.548 0.585 242 89 54 
23.6 1.39 Loilg 402 191 144 
25 0.94 1.0 496 214 69 
Dy 1.96 1.78 395 158 121 
30.4 0.84 0.77 291 95 84 
33.8 0.736 0.743 328 58 43 
a 478 195 110 
U 292 84 45 
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lengths (5 460-5 470 A), showing good agreement. The 
differences in the values of » may be caused by several 
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factors. The uncertainty in the J-determination (the 
thickness of the crystals has been supposed to be equal 
to the horizontal diameter) ; neglecting the change of 
refractive index with the wavelength may also cause 
an error ; the interference effects, too, are neglected ; 
they are dependent on the wavelength. The most 
important factor is the real path of light within the 
crystal. This problem was mentioned also by |AntTo- 
Nov-Romanovsky [3] who distinguished the per- 
fectly regular shaped crystals from the irregular ones. 
Through a perfectly regular-shaped crystal, ‘e.g., a 
sphere or a cube, the light passes without internal total 
reflection, while in an irregular shaped crystal rays 
passing through are repeatedly totally reflected. The 
true value of the absorption coefficient is not the ave- 
rage value calculated from the measurements, u, but 
their minimum value. It may be assumed that the 
minimum values are those found for crystals through 
which the radiation passes only once. This corresponds 
to the perfectly regular-shaped particles. The true 
value of uw cannot be furnished by any statistical 
method. Single microcrystals examined for brightness 
justify this supposition. 

The microscope photometer and the empirical model 
of evaluation were experimentally checked on K,Cr,O, 
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Fic. 3. — The light emission of microcrystals of different particle size. 
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crystals. A very good agreement was found for the 
5 460 A absorption coefficient determined for micro- 
crystals and large crystals. 


Relation between light emission and absorption of 
luminescent microcrystals. — In 1950, Bopo [14 ] exa- 
mined the surface brightness of halophosphate micro- 
crystals under 2537 A excitation by a visual method. 
He found arelation between the average surface bright- 
ness and the ultraviolet absorption of the crystals. 

A slight modification of the experimental arrange- 
ment of figure 1 makes possible the measurement of the 
light emission of single crystals. Figure 3 shows the 
results of measurements with the 90 ZnS-10 CdS-Cu 
sample. The light emission increases with increasing 
particle size, but no continuous curve has been found ; 
the experimental points are situated in an area limited 
by two curves. The upper limiting curve may be des- 
cribed by the equation : 


Decca keh® 
k being a constant. The lower curve is : 
I =k(P/l —e) 


In this equation p is equal to the minimum value of 
the 3 650 A absorption coefficient, determined by the 
microscope-photometer. The upper limiting curve corres- 
ponds to the irregular shaped crystals. The light emis- 
sion is proportional to the absorbed ultraviolet energy. 
Ultraviolet radiation is practically perfectly absorbed 
in consequence of multiple internal total reflections. 
The absorbed energy is proportional to the crystal 
surface, /?. The lower limiting curve on figure 3 repre- 
sents the case of perfectly regular shaped crystals. In 
practice, the majority of the crystals fall between the 
two limiting cases. Using » in the formula gives the 
dotted line in figure 3, corresponding to the average 
absorption. 


Experimental results. — Table II contains the 
absorption coefficients for some ZnS-type materials. 
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The crystals examined were prepared by I. Hanaos 
and Mrs K. Toru. 


TABLE II 
3650 A 4060 A 4358 A 
Specimen = cy a3 

U U U. U U. ue 
hex. ZnS-Ag sample 1 ...| 372]120 
hex. ZnS-Ag sample 2. ...}| 538/148 
hex. ZnS-Ag sample 3. ...} 397/124 
hex. ZnS sample 3 

WithouteNcken vente 340) 57.3 

Cuba Zn SoNe are eerenvereds 364|139 
laen ss PARIS) Cll eso no ofa. coon ADS 220 VAAN S77 |e Ooms 
hex. 90 ZnS-10CdS-Cu...| 478}292 195} 84] 110) 45 
hex. 50 ZnS-50CdS-Ag ...| 843/636 740} 555] 131} 80 


Reflectivity of microcrystals. — A modification of 
the optical arrangement of figure 1 (vertica] illumina- 
tion) made reflectivity measurements on single micro- 
erystals possible. From the reflectivity data, the refrac- 
tive index could be estimated. 

A full description of the apparatus, the experimental 
methods, the optical properties of spherical micro- 
crystal and detailed data on ZnS-type microcrystals 
will be given in the Acta Phys. Hung. 


Acknowledgement. — The author wishes to express 
his thanks to his co-workers, J. Anam, I. Haneos and 
Mrs K. Torn. 
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The dependence of the human electroretinogram on the shape 
of the stimulus as a function of time * 


Lucia Roncur and Silvano Graz 
Istituto Nazionale di Ottica, Arcetri, Firenze (Italie). 


Summary. — The human electroretinogram has been studied for various time-gradients of illumination of the retina, using an intense 
white stimulus, for both a light and dark adapted eye. In general the cone response appears to be masked by rod activity. By 
varying the time of rise of the stimulus between 1 and 70 msec no remarkable variations in the ERG are evident. On the other hand, 
when the time of rise ranges from 70 to 600 msec remarkable changes in both shape and size occur, and the latency time increases. 


SOMMAIRE, — L’électrorétino gramme (ERG) a été étudié pour différents gradients temporels d’éclairement sur la rétine de ’ eil humain, 
éclairé par une lumieére blanche intense, adapté a V’obscurité ou non. La variation du temps d’établissement du stimulus entre 1 et 


70 msec, n’a pas apporté de changement notable dans VERG. 


D’autre part, quand ce temps passe de 70 & 600 msec, des changements remarquables se produisent dans la forme et la grandeur 


de VERG, et la latence augmente. 


ZUSAMMENFASSUNG. — Das Elecktro-Retinogramm wurde fiir das menschliche Auge bei Hell- und Dunkeladaptation untersucht 
wenn bei einem starken weissen Reiz verschieden schnelle zeitliche Anderungen der Beleuchtungsstdrke auf der Netzhaut benutzt 


(*) This research has been made possible through the support and sponsorship of the Air Research and Development Com- 
mand, United States Air Force, between its European Office and the Istituto Nazionale di Ottica under content N° AF61 


(514)-634 C. 
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werden. Bei zeitlichen Anderungen der Reizzunahme zwischen 1 und 70 msec. ist kein Unterschied in den Elektro-retinogramm 
zu bemerken. Wenn aber die Zeit fiir das Anwachsen des Reizes zwischen 70 und 600 msec. liegt, so ergeben sich betrdtliche 
Abweichungen in Form und Grésse und die Latenzzeit wéchst an. 


Introduction. — The electroretinogram (ERG), 
represents the only measurable objective index of 
what takes place in the retina in response to change 
in illumination [4]. The “ changes in illumination ” 
are usually realized by means of a light stimulus the 
form of which is a very complicated function of time. 
In a simple case, the curve presents four angular 
points, A, B, C, D, as represented in figure 1. 


A D 


Fig, 1. — A simple light stimulus (intensity 
as a function of time). 


It is conceivable that if the duration of the stimulus 
is sufficient, we may obtain two well-separated res- 
ponses both at the beginning and at the end of the 
light stimulus (the so called on-effect and off-effect 
respectively). Also it is conceivable that the on-effect 
is a response depending on the variation A only, pro- 
vided the slope of the segment AB is suitable. 

These problems involve the study of two factors 
influencing the ERG: the total duration of illumina- 
tion and the shape of the light stimulus. The former 
was taken into account for the first time by pe Haas 
in 1903 [2], and until Wirarx [3] few authors have 
attempted a similar investigation. On the other hand, 
the shape of the stimulus as a factor has never been 
examined, so far as we know. 

We have undertaken a program of research begin- 
ning with very simple light stimuli, containing only 
the first two variations A, B represented in figure 1, 
i.e., we used time gradients of illumination (fig. 2) 
and our aim was to examine how the ERG is influen- 
ced by varying the slope of the gradient. 

Our research originated at the suggestion of Pro- 
fessor Toratpo who pointed out similarities in the 


A 


Fic. 2. — Time gradient of illumination 
(intensity as a function of time). 


Illumination in metre candles 


0.05 0.10 0.15 0.20 0.25 0.30 
Time in seconds 
Fic. 3. — The development of the visual response 
during stimulation with white light (Broca and 
SULZER). 


results of widely different experiments, observed by 
some authors. For instance, in 1945, Wriaut [ 4] said 
that “‘ the general similarity between the curves repre- 
senting the Broca-SuLzerR phenomenon (fig. 3) and 
the early part of the ERG should be noted ”’. On the 
other hand, ToraLpo (personal communication) notes 
the similarity between the early part of the ERG and 
the behaviour of the differential threshold in a field 
with a spatial gradient of illumination [5 ] (see fig. 4). 
The connecting link between the spatial and temporal 
phenomena is assumed to be in the small involuntary 
movements of the eye. 


A 


threshold 


Differential 


Visual 


angle 


Fic. 4. — The differential threshold (full line) mea- 
sured on a field with a gradient of luminance (broken 
line); the abscissae are the visual angles measured 
from the centre point of the field (F1oRENTINI). 
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Fic. 5. — Stimulating system (for significance of symbols see text). 


The crucial point is the comparison between the 
subjective response and the objective response (ERG). 
A more specific connection might be derived by consi- 
dering the relation between the b-wave and the rise 
in frequency of the optic nerve discharge. (Such an 
analysis would necessarily assume a similarity in 
behaviour of the nerve fibres of cold- and warm- 
blooded animals). This rise in frequency is similar to 
the rise in the sensation of brightness in man (BRocaA 
and SuLzeErR effect) [6]. We hope that the research 
which we have begun will be a contribution to the 
question ‘‘ how is the form of the stimulus connected 
with the form of the ERG ? ”; in other words, our 
purpose consists in clarifying under what conditions 
the eye gives an electrical response which is indepen- 
dent of the energy-time distribution of the light sti- 
mulus, and under what conditions, on the other hand, 
the ERG is influenced by the energy-time distribu- 
tion. 


1. Apparatus. The apparatus used in the study 
of the electroretinograms consists of two main parts ; 
one of these, a stimulating system, is used to deliver 
time gradients of illumination to the eye and the other, 
a recording system, is used to trace the electrical res- 
ponse. 


Stimulating system. — Figure 5 provides a view of 
the principal elements of this system; an objective 
Ob forms at D the image of a 6 amp., 6 volts ribbon 
filament lamp, and the lens L, placed close to the eye, 
projects the image to infinity, and images the aperture 
of Ob at the eye pupil. 


Faraday 


Cage 


The observer’s eye views the filament through L ; 
the diameter of the light beam entering the eye pupil 
does not exceed three millimeters. Great care is there- 
fore required in order to assure immobility of the 
observer’s position. Fixing the head is by the usual 
dental impression method. The observer is further 
helped in avoiding voluntary eye movements by a 
dim red fixation point placed at the centre of the 
screen D. 

The stimulus has the form shown in figure 2 ; as can 
be seen, the exposure time is not limited. A slotted 
plate guided by vertical rails falls under gravity from 
a given height and allows a linear time gradient of 
illumination to be realised. The time of rise ¢, of the 
light stimulus depends on the drop-height ; the rising 
of illumination is linear if the aperture of Ob has a 
particular form, which depends also on the drop- 
height. 

This method has been adopted for ¢t) ranging from 
1 to 60 msec. For obtaining values of ft) greater than 
60 msec., the aperture of Ob is made rectangular and 
its size is independent of ¢, ; the falling plate is driven 
by a synchronous motor, whose speed may be varied. 


Recording system. — A block diagram of the recor- 
ding system is shown in figure 6. Electrical potentials 
are recorded between an active electrode mounted in 
the corneal bulge of a contact lens and an indifferent 
electrode worn on the forehead. Signals picked up by 
these electrodes are taken to a capacitance-coupled 
amplifier, the output of which is connected to a cathode 
ray oscilloscope. Permanent records are made by pho- 


Camera 


Pia. 6. — Block diagram of the recording system. 
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the full line refers to light adaptation, the broken line 
of the stimulus is indicated in each figure. Right: the form of the light stimuli. 
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tographing the trace on the oscilloscope. The gain of 
the amplifier may be varied from 600 to 1 000. 

The sweep is provided externally ; it consists of a 
microswitch placed conveniently on the vertical rails ; 
the falling plate acting on the microswitch provides 
the current necessary for the displacement of the spot 
along the time axis. 

Besides the channel providing a complete record of 
all activity appearing between the electrodes, another 
channel has been used in order to provide a complete 
record of the form of the light stimulus. This channel 
is also capacitance coupled and is fed from a photocell 
which receives a part of the light stimulus by reflection 
from a mirror placed on the screen D shown in figure 5. 


2. Procedure. — During an experimental session, 
the subject sits in a FaArapaAy cage and first dark or 
light-adapts for fifteen minutes. When dark-adapta- 
tion is required, no light passes through the walls of 
the screened cage, but when light adaptation is requi- 
red, the walls are diffusely illuminated at a given level. 

A few seconds before a test-stimulus is to be admi- 
nistered, the experimenter signals the observer, who 
moves his eye to the central fixation point, the other 
eye being covered. After the stimulation, the eye is 
again light or dark-adapted as required. The stimuli 
are administered at five minute intervals. 

The angular aperture of the beam entering the eye 
pupil is nearly one degree. The illumination at the eye 
is 100 lux during stimulation, 15 lux during light- 
adaptation. Two normal observers took part in these 
experiments. The contact lenses were produced by 
Firm ‘“‘ Opticon ” in Milan. 


3. Results. — The experimental results given by 
observer S. G. are represented in figures 7 to 11. The 
total number of responses examined was 834. Figure 7 
shows the average of the responses obtained for the 
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Fic. 8. — The height of the b-waye, in wu V, is plotted against ty, 
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times of rise of the stimulus, ¢), indicated for each 
curve. The broken line refers to dark-adaptation 
(d.a.), the full line to light adaptation (l.a.) corres- 
ponding to an illumination at the eye of 15 lux. The 
form of each light stimulus is represented at the right 
of the corresponding response. 

From the figures we see that when the rate of rise 
of the stimulus decreases, the rate of rise of the res- 
ponse (b-wave) decreases also. After the peak, when 
ty) is large, there is no return to the base line, and we 
could suggest, perhaps, the presence of another com- 
ponent of the b-wave. 

In addition, a curious effect has been noticed ; for a 
given value of ¢) greater than 150 msec., we obtained 
many unfamiliar responses, consisting for instance of 
two or more subsequent b-waves. In the last curve of 
figure 7 we see for instance that, in light adaptation, 
the response may consist either of a slow rise of poten- 
tial, or in a line parallel to the time axis. 

In figure 8 the height of the b-wave is plotted against 
to. In figure 9 the rate of rise s,, of the b-wave is plotted 
against t¢), and figure 10 represents the behaviour 
of s, as a function of log 4/t. It is evident that 1/t 
represents the rate of rise of the light stimulus, if / 
expresses the luminous level in arbitrary units. In 
figure 11 the latency time of the b-wave is plotted 
against fo. 

Let us consider now the other observer L. R., from 
whom about 100 records have been obtained in the 
light adapted state. These results seem in general 
agreement with those of S. G. However we may note 
the following differences : the heights of the b-wave 
are nearly 40% greater, and the latency times nearly 
50°% less than in the records of the first observer. 
Besides, for L. R., a number of small waves are detec- 
ted superimposed on the b-wave. Figure 12 shows 
three records each referred to the value of tj indicated 
on the right. 
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the time of rise of the stimulus, expressed in seconds. 


bserver S. G. (d. a., dark adaptation ; 1. a., light adaptation). 
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4. Discussion. — We may attribute the above men: 
tioned differences in the form of the response to a 
given light stimulus to disturbances such as fatigue 
due to prolonged fixation. However, we read in GRa- 
nit [7]: ‘‘ ... Sometimes much to the experimenter’s 
surprise, two quite different responses to the same sti- 
mulus may alternate, or even appear or disappear, 
according to no obvious rule... The eye seems to 
remain in excellent condition during the experiment 
and shows no sign of fatigue ”. In his opinion, the 
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Fic. 10. — The decimal logarithm of the slope of 
the b-wave is plotted against the logarithm of the slope 
of the light stimulus. 


erratic behaviour of the component PII (responsible 
for the b-wave) must account for a phenomenon simi- 
lar to the rotation activity among the elements parti- 
cipating in the nerve discharge. In other words, these 
experiments suggest that the retinal elements which 
respond to light with the reaction PII do not all 
behave alike ; PII therefore shows a certain degree of 
complexity. Although these effects are very frequent 
at low intensities, they cannot be excluded at high 
intensities. Analogously we consider that in our case 
the above effect is very frequent for large values of fo, 
but sometimes it has been observed also for small 
values of ft. 

By examining our records, we notice that for ¢) less 
than 70 msec., the b-wave starts before the ‘“ plateau” 
in the light stimulus is reached. On the other hand, 
for t) greater than 70 msec., the response seems to be 
due to both variations of illumination A and B (fig. 2) 
occurring in the stimulus. For times of rise of the sti- 
mulus ranging between 1 and 70 msec., the response 
of the eye seems not to be influenced by variations 
of t; on the other hand, when ¢, is greater than 
70 msec., the height and the rate of rise of the b-wave 
decrease, and the latency time increases; but the 
decrease of the rate of rise is more pronounced than 
that of the height. 

Let us consider now figure 10; the curve obtained 
in dark adaptation may be divided into two main 
parts ; one of these is independent of t, ; the other one 
may be considered as a straight line, the slope of 
which is unity. This means that the rate of rise of the 
b-wave is a linear function of the rate of rise of the 
light stimulus. The behaviour of the curve referred to 
light adaptation seems to be more complicated. 

Further, we may add that in the dark adapted state, 
component PII is dominant with respect to light adap- 
tation. It seems reasonable to ask what has happened 
to the cone response in the records given by observer 
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Fig. 11. — The latency time of the b-wave is plotted againts the time of rise of the stimulus. 


S. G. For é) ranging from | to 10 msec, the a-wave is 
very small (10 to 20 V) and for ¢ greater than 10 msec. 
the a-wave is not detected. 

For ? less than LO msec., a very small wave is some- 
times indicated, suggesting a reduced form of Moro- 
KAWA and Mita’s X-wave [8]. But as the curve is of 
the same order of magnitude as the noise level, no 
definite statement can be made at this stage. 

The intensity of our stimuli may be considered high, 
and in light-adaptation the cone response should be 


Obs: L.R 
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100 msec 


Fig. 12. — Some records obtained by observer L. R. 


uncovered *, according to Granit [9] [10]. The 
suppression of cone activity by the rods might depend 
on the length of time of illumination. According to 
Wirth [3], the cone response in light adaptation 
when strong stimuli are used is elicited when the expo- 
sure time is less than 20 msec. The same limit holds 
also in dark-adaptation. In our case, by using time 
gradients of illumination, the exposure time is not 
limited ; the lack of cone activity in the responses 
given by S. G. therefore agrees with Wirtn’s results. 

In the responses given by L. R., the a-wave is very 
small, but the cone response does not seem to be com- 
pletely suppressed. A small wave near the beginning 
of the b-wave is present in all the responses obtained 
for fy ranging between {1 and 300 msec. The X-wave 
seems to appears for fy = 10 msee., while, for 4 = 1 
msec., the b-wave has a composite nature. The 
behaviour of L. R. is not in disagreement with 
Wirtn’s results. We suggest that if the light stimulus 
is very short and rises rapidly, then only the cone 
system is able to react. 

If the light stimulus is long and rises slowly (within 
limits), the rods response may be recorded. If the sti- 
mulus is rising rapidly and of long duration, the res- 
ponses of both systems may coexist, or as in the case 
of S. G., the rod response being dominant suppresses 
cone activity. It should depend on the relation bet- 
ween the time-constants of both cone and rod systems 
in the retina of each observer. 


Conclusion. In the present study we examined 
the influence on the human electroretinogram of 
intense light stimuli, having time gradients of illumi- 
nation. 

By varying the rate of rise of the stimulus, both the 
shape and size of the electroretinogram undergo modi- 
fications if the time of rise is greater than 100 msec. 
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When the stimulus rises slowly, the response is of 


the low-intensity type. When it rises rapidly, the res- 


ponses of one observer probably show signs of cone a 
activity ; on the other hand, in the corresponding [2] 
responses of the other observer, the cone activity 16] 
seems to be completely masked by that of the rods. ye 
[5] 
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The effect of defocusing and third order spherical aberration 
on the transfer function of a two dimensional optical system 


G. B. Parrent and C. J. DRANE 
Air Force Cambridge Research Center. Bedford, Mass., U.S.A. 


The effect of defocusing and third order spherical 
aberration on the transfer function of a two-dimen- 
sional optical system is the subject of a study now 
reaching completion at this Laboratory. Some results 
of general interest are already apparent and form the 
subject of this communication. 

The recent emphasis on the optical transfer func- 
tion as an analytic tool precludes the need for presen- 
ting an extensive mathematical development. It is 
readily shown that for incoherent illumination the 
transfer function <(u, v) is the convolution of the aper- 
ture illumination t(y, v) with its complex conjugate 


(Duffieux) [ 4 ]. 


(a) 


Prt ge 


(1) (u,v) = l| t(u', v') t(p.’—p,v'—v) du’ dv’. 


ew 


where (uy, v) are the aperture variables. 
With uniform transmission over the aperture 


a, Vv) = eiA(y, v) 
= (0) 


(2) v,¥ C aperture 


vu, ¥ © aperture 


where A describes the deviation of the actual wave- 
front from a reference sphere. For third order spheri- 
cal aberration and defocusing 


(3) 


where a and b are the coefficients of spherical aberra- 
tion and defocusing, respectively. After a suitable 
change of variables equation (1) becomes 
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— Comparison of theoretical and experimental transfer functions for one wavelength 


of spherical aberration in two dimensions. 


b — 2 focal shift 
a total longitudinal aberration 
= focal setting. 


where ¢ = 


The inner integral will be recognized as Arry’s In- 
complete Integral. However, existing tables were 
inadequate and the integrals were evaluated numeri- 
cally. Simple generative relationships are readily 
obtained which make possible the extension to larger 
aberrations. 

Since the convolution integral represents a smoothing 
action, it might be expected that the chief difference 
between the one and the two-dimensional function 
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would be the extra smoothing effect in the two-dimen- 
sional case, that is, changes in the one-dimensional 
functions (PARRENT) [2] would be more pronounced 
than corresponding changes in the two-dimensional 
functions. That this is indeed the case is seen from 
figure 1. The excellent quantitative agreement between 
our calculated results and those obtained on the Ana- 
log Interferometer at Boston University* is shown in 
figure 2. 

As would be expected from geometric optics, a much 
better system is indicated by the midfocus curves in 
figures 1 and 2 than by the paraxial or marginal focus 
curves, since the extremes show reduced contrast and 
considerable amounts of spurious resolution. A detai- 
led analysis, however, demonstrates that no single 
focal setting is best for all spatial frequencies. If, the- 
refore, the best focus for periodic detail were to be 
defined as that focal setting at. which the transfer 
function is maximum for the desired frequency, it 
would be possible to plot best focus versus frequency 
as in figure 3. As would be expected from the « smoo- 
thing » integral, the shift in the two-dimensional case 
is less pronounced than that observed in the one dimen- 
sional case. This curve is in good qualitative agree- 
ment with the experimental results reported indepen- 
dently by Kineastake [3], MacDonaip [4], and 
Scuack and Batt [5]. 

In most applications of optical systems, one observes 
random distributions, not sine wave targets. A perti- 
nent question therefore is where lies the best focus 
for random detail. Since the root mean square fluc- 
tuation in the image is related to the square of the 
transfer function, a reasonable criterion (in the absence 
of spatial phase shift) is the area under the squared 
transfer function (Schade’s equivalent bandpass) [6] : 


* A device which simulates a shearing interferometer by 
making use of Morre fringes. 
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between paraxial (0) and marginal (1) for one wave- 
lenght of spherical aberration in one and two dimen- 
sions. 
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Since spherical aberration does not involve a spatial 
phase shift, it is reasonable to define best focus for 
random detail as that focal setting for which the band- 
pass is a maximum. In figure 4 the equivalent bandpass 
versus focal setting for one wavelength of spherical 
aberration is plotted for a one dimensional and a 
two-dimensional optical system. It should be noted 
that the maxima of these curves are shifted toward 
the paraxial focus, which is in qualitative agreement 
with the results of Lansravux [7 ]. Further, it is signi- 
ficant that the shift is more pronounced for the one- 
dimensional than for the two-dimensional system. 
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The calorimetric determination of the diffuse reflectance of magnesium oxide. 


Gy. GERGELY, 
Research Institute for Telecommunication, Budapest. 


The diffuse reflectance of magnesium oxide has been 
studied by many authors. Recently Benrorp, Knowles 
MippLeton, Scuwarz, Lioyp [1-4], TeL_ex and 
Watpron [5] described the results of highly precise 
measurements. All of them used some modification of 
the integrating sphere method. The results of the 
different measurements agreed well. Here we do not 
wish to mention more than two aspects of the problem. 
The above mentioned authors found that the diffuse 
reflection of magnesium oxide is affected by the con- 
ditions of preparing the coatings, further it varies 
with irradiation, age, etc. of the coating — as descri- 
bed by Jacquez and co-workers [6]. The changes of 
the reflectance especially in the ultraviolet part of the 
spectrum cannot be neglected. Last year TELLEX, 
Watpron and Hammonp [7] drew attention to the 
importance of the thickness of the coating. 

In order to reduce the errors in the determination 
of the reflection spectra of diffuse powders, e. g., lumi- 
nescent powders, a comparison standard was used in 
our laboratory. Our standard diffuse reflectance sur- 


face was a layer of magnesium oxide powder 1 mm 
thick, in order to avoid the uncertainties of the magne- 
sium oxide coating. An 1 mm thick layer is an opti- 
cally infinite thick layer. It can easily be prepared and 
changed. 


The preparation of the magnesium oxide layer. — 
The preparation of thick magnesium oxide coatings, 
especially the coating of an integrating sphere having 
thick layer, can be a very difficult task. TELLEXx and 
Watpron described a procedure, but their method 
is a difficult one. A thick magnesium oxide layer can 
be easily produced by spraying magnesium oxide 
powder on a sample holder plate. We used pure magne- 
sium oxide powder for this purpose. The powder was 
dehydrated at 1 000°C as described by R. Razouk 
and Mrxualit [ 12]. 


The calorimetric determination of the diffuse reflec- 
tion coefficient. — The diffuse reflectance of thick 
magnesium oxide powder layers cannot be determined 
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Fic. 1. — Experimental arrangement. 1 : 
4: magnesium oxide sample ; 5: blackened plate ; 6: 
(ScHoTtT) ; 9: 


by the integrating sphere method used by the above 
quoted authors [1-5]. The reflectance of a non-lumi- 
nescent powder can be easily determined by a calori- 
metric method. ALENcEv [8] and Bopo [9] used a 
microcalorimeter for the determination of lumines- 
cence efficiency. We improved ALENCEv’s and Bop0o’s 
experimental technique. Instead of thermocouples, 
thermistors were used, as it was done by Hutcuinson 
and WuirteE [ 10]. 

Figure 1 shows the sketch of our apparatus. The 
magnesium oxide sample to be examined was placed 
on a sample holder plate, within an ultrathermostat 
of the HorppLer type. Following Bop5’s method, 
another sample holder plate, blackened with carbon 
black, was used for the determination of the intensity 
of radiation. Both sample holder plates were placed 
upon the copper plate of a disc thermistor connected 
with another thermistor by a thin copper rod. The 
two thermistors formed the arms of a bridge. 

The magnesium oxide sample and the carbon black 
surface were irradiated by a high pressure mercury 
lamp. The strong spectral lines of the lamp were selec- 
ted by Scuort interference filters. The radiation entered 
the thermostat through a quartz absorption cell, 
through which the distilled water of the thermostat 
was circulated. 

Irradiating the magnesium oxide and the carbon 
black samples raised their temperature. The tempera- 
ture rise was detected by the thermistors by means of an 


ultra-thermostat ; 2 : thermal insulator frame ; 3 


quartz absorption cell; 7 : condenser lens ; 8 
high pressure mercury lamp ; 10: electricon thermistor bridge ; 11 : electronic voltmeter. 


: disc thermistors ; 
: interference filter 


electronic bridge developed by Gy. ArmAssy. This 
apparatus will be described in another paper [ 14 ]. The 
method of measurement was similar to that of Bono, 
used for luminescence efficiency determination. The 
formulae derived by Bond will be employed also in 
this paper. 

Irradiating the magnesium oxide sample, its tem- 
perature rise is given by the differential equation : 


(1) KdT = I9,(4 —rmy) dt — B,(T — T,)at . 
where 

T is the temperature of the sample, 

T, is the ambient temperature, 


K is the heat capacity of the system, 

t is the time, 

Z,,1s the intensity of radiation at the wavelength 4, 
r, 18 the reflectance of MgO, 

B, is the heat transfer coefficient. 


Solving equation (1) gives 
vf it here ») B 
(oy erie pee a (1 Sep zt): 


A similar equation gives the temperature rise of the 
carbon black sample. Using a similar notation, where 
ris the reflectance of the carbon black, gives 


J deatin) B 
ia a By E exp— Set]. 
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Reaching the thermal equilibrium at approxima- 
tely t = 10 minutes : 


al 


ekp += he expe 0: 
Thus, we get 
Le 1 (1 =) tee ’ 
2&2, 
(4) 
L(A Ley 15%, 
M,= ‘Fp > Cie = We ) 


denoting by M,, and M, the assymptotic value of 
temperature rise, by 1/<, and 1/<, the time constants 
of the system, irradiating the MgO, and the carbon 
black respectively. M,, M,, c¢, and ¢, were determined 
by measuring the rate of heating of the system. It was 
found that r, ~ 0,015 independently of the wave- 


length. 


Results. —- Table | contains the result for some 
intense mercury lines. The 2537 A radiation was pro- 
duced by means of a low pressure mercury discharge 
lamp. For the purpose of comparison, Table I con- 
tains the results of Benrorp and co-workers. It may 
be seen that agreement was found with the result 
published in [ 1-6 }. In the ultraviolet part of the spec- 
trum, the reflectance of our sample was less than that 
found by Benrorp and co-workers. This difference 
may be explained by the small traces of iron impurity 
of our MgO powder. 
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thermistor bridge, to Mr Gy. ALmAssy, for the help in 
the measurements, to Mr J. Anim, for building the 
electronic device, to Mr P. Nagy, and, for building the 


microcalorimeter, to Mr E. NéMetu. 


TABLE | 
Diffuse reflectance of magnesium oxide 
Wavelength 

calorimetric BENnForp’s and 

determination co-workers’ results 
2 537 oo % 945% 
3 650 91.2% 94.5% 
4 060 SW YA Mel ah 
4 358 96.4 % 98 % 
5 460 99 % 99 % 
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Collogue international sur les Problémes physiques de la télévision en couleurs 
(Paris, juillet 1957) 


Sous le patronage de |’Union Internationale de 
Physique Pure et Appliquée, un Colloque consacré a 
Etude des Problémes Physiques que pose la Télévi- 
sion en Couleurs sera organisé 4 Paris du 2 au 6 juillet 
1957. Les séances de travail auront lieu au Conserva- 
toire National des Arts et Métiers, 292 rue Saint- 
Martin a Paris. 

Programme provisoire. — Au cours de P Assemblée 
pléniére du Comité Consultatif International des Radio- 
communications tenue a Varsovie en aout-septembre 
1956, la Commission XI (Télévision) a recommandé 
la recherche de méthodes et de résultats quantitatifs 
relatifs a l’évaluation des qualités des images de Télé- 
vision et approuvé la constitution d’un groupe de travail 
a cette fin. Les physiciens, les ingénieurs électroniciens 


et opticiens, les spécialistes de la télévision trouveront 
ainsi ’occasion de se réunir et d’étudier les points 
suivants : 

A. Comportement de lceil, 

B. Prise de vues et restitution de l'image, 

C. Mesure des résultats obtenus dans la restitution 
d’une image colorée, 

D. Systeme de codage utilisés pour la transmission 
des signaux de Télévision en Couleurs, : 

Communications. — Les communications a présen- 
ter au Colloque devront étre annoncées trés prochai- 
nement et en tous cas avant le 1° mai 1957 pour 
permettre l’impression et la distribution d’une liste 
de titres et de résumés en temps utile. 
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Les communications pourront ¢tre présentées et 
commentées dans la langue choisie par lauteur ; 
toutefois Pusage du franeais et de Vanglais sont 
recommandes. 


Démonstrations, visites, excursions. — Enfin des 
démonstrations, visites et excusions seront organisées. 
La cotisation du Colloque sera de lordre de 
| 500 Francs Franc¢ais pour chaque participant. Etant 


INFORMATIONS 


[o pr. Acta 


donné Vaffluence des touristes & Paris en juillet, le 
Comité d’organisation sera reconnaissant a tous ceux 
qui envisagent de participer au Colloque de le faire 
savoir dés que possible a l’adresse ci-dessous : 


Colloque International sur les Problemes Physiques 
de la Télévision en Couleurs, Conservatoire National 
des Arts & Métiers, 292, rue Saint-Martin, Paris, 3°. 
France. 


ERRATA 


« The defocused image of sinusoidal gratings » by W. H. Steel, Optica Acta, 2, 1956, 65. 


Page 68 : 
Right hand line below Table 1 
fOr aOsORe 


Ordinate scales, Figure 2 


3 
for s(2ecoe" aes: 


read "—0.5". 


3 
‘read “(2.00128 one 


for "0.05", read “—0.05 ". 


Page 69: 
Top line, equation (14) tor 225) | COS, 7a) pened | te | COs ae 
(63) (e3) 
Pages 70: 
2 i 2h eee 2 
Penultimate expression for" “_— ", read " pa Mee i Me 
2 wo 2 wo 
Pages 73 : 
Second line of Appendix for " « The frequency of a... ". 
read " « The frequency response of a... ”. 
Pages 74: 
Reference [12] for " 1954" read ” 1953". 
Throughout (pp. 65, 66, 68, 
69, 70 and 73) for " spacial” read " spatial ’. 


ADDENDUM 


J. Fockxr. Wellenoptische Untersuchungen zum Offnungsfehlers, Optica Acta, 3, 1956, 126. 


Literaturhinweise : 


Siehe auch Prcur, J., Optische Abbildung, Braunschweig, 1931. 
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